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ABSTRACT 


This  report  covers  the  data  reduction  phase  of  the  project. 

The  work  was  concerned  with  the  determination  of  the  payload  attitude  of  a  number 
of  hlgh-altltude  rocket  flights.  From  this  Information  the  pointing  direction  of  the  vari¬ 
ous  directional  detectors  aboard  the  rockets  is  available  for  analysis.  The  telemetry 
data,  transmitted  from  the  rookets,  was  prepared  in  a  digital  form  on  tape  for  input  to 
a  computer  program.  The  computer  program  calculated  the  payload  attitude  in  a  John¬ 
ston  Island  coordinate  system.  The  outputs  of  the  detector  functional  values  and  detector 
orientation  parameters,  trajectory  parameters,  and  associated  telemetry  data  are  avail¬ 
able  in  printed  form  as  a  function  of  time  In  Increments  of  11  milliseconds. 

The  report  discusses  the  theoretical  basis  for  attitude  determination  of  the  payloads 
and  presents  the  computer  program  in  its  final  form  as  used  In  the  data  reduction.  A 
description  of  the  data  reduction  results  is  Included  as  the  laBt  section  and  discusses 
the  final  form  of  the  computer  outputs  and  character  of  the  data. 

Binary  coded  digital  data  tapes  were  prepared  during  the  computer  processing  of 
the  data  for  use  as  Input  to  a  data  analysis  program  currently  in  progress. 

The  results  and  analysis  of  the  Project  6.2  Fish  Bowl  data  will  appear  In  the  Final 
Report  of  Contract  DA-49-146-XZ-201. 
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CHAPTER  1 
INTRODUCTION 

This  report  covers  the  data  reduction  phase  of  Project  6.2,  and  is  Intended  as 
an  addendum  to  POR-2017  (Reference  1). 

The  work  was  performed  under  DASA  Contract  DA-49-1 46-XZ-123  (EOS  ln- 
house  work  authorization  WA  2194).  It  was  concerned  with  the  determination  of 
the  payload  attitude  of  a  number  of  high-altitude  rocket  flights.  From  this  Infor¬ 
mation  the  pointing  direction  of  the  various  directional  detectors  aboard  the 
rockets  Is  available  for  analysis.  The  telemetry  data,  transmitted  from  the 
rockets,  was  prepared  In  a  digital  form  on  tape  for  Input  to  a  computer  program. 
The  computer  program  calculated  the  payload  attitude  In  a  Johnston  Island  coor¬ 
dinate  system.  The  outputs  of  the  detector  functional  values  and  detector  orienta¬ 
tion  parameters,  trajectory  parameters,  and  associated  telemetry  data  are 
available  In  printed  form  as  a  function  of  time  In  increments  of  11  milliseconds. 

The  report  discusses  the  theoretical  basis  for  attitude  determination  of  the 
payloads  and  presents  the  computer  program  In  its  final  form  as  used  In  the  data 
reduction.  A  description  of  the  data  reduction  results  is  Included  as  the  last  sec¬ 
tion  and  discusses  the  final  form  of  the  computer  outputs  and  character  of  the  data. 

Binary  coded  digital  data  tapes  were  prepared  during  the  computer  process¬ 
ing  of  the  data  for  use  as  Input  to  a  data  analysis  program  (Contract  DA-49-1 46- 
XZ-201)  currently  In  progress. 

The  presentation  of  the  results  and  analysis  of  the  Project  6.2  Fish  Bowl  data 
will  appear  In  the  Final  Report  of  Contract  DA-49-146-XZ-201. 
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CHAPTER  2 

THEORETICAL  DISCUSSION  OF  ATTITUDE  DETERMINATION 
2.1  DEFINITION  OF  PROBLEM 

To  adequately  describe  the  debris  cloud  resulting  from  the  nuclear 
explosion,  it  is  essential  to  know  the  direction,  at  any  time,  of  each 
detector  on  the  payload  relative  to  an  inertial  coordinate  system. 

Since  the  payload  flight  time  is  small  relative  to  the  earth's  period 
of  rotation,  the  coordinate  system  can,  for  most  practical  purposes, 
be  fixed  in  the  earth  without  resulting  in  a  significant  error.  Thus, 
during  an  interval  of  400  seconds,  which  is  the  average  flight  time, 
the  earth  will  rotate  through  an  angle  of  1.68  degrees.  This  error  is 
probably  less  than  the  accuracy  with  which  one  can  determine  the  pay- 
load  attitude  by  means  of  the  magnetometer  data.  A  convenient  coordi¬ 
nate  system  for  determining  the  payload  attitude  is  the  x,  y,  z  system 
shown  in  Figure  2.1(a)  which  has  its  origin  at  Johnston  Island.  The  z-axis 
points  along  the  local  vertical  (out  of  the  earth),  the  y-axis  points 
to  true  (geographic)  north,  and  the  x-axis  points  due  east.  Here,  D 
denotes  the  sensitive  axis  of  a  detector  (e.g.,  a  gamma  scanner)  where 
9  and  4  are  the  polar  coordinates  of  this  axis.  Figure  2.1(b)  defines 
the  fixed  polar  angles  0',  4’  of  the  detector  axis  relative  to  the  x1 , 
y’  and  z’  magnetometer  axes.  The  values  of  0',  4'  for  the  four  detec¬ 
tors  in  the  payload  are  shown  in  Table  2.1.  Although  the  origin  of  the 
x',  y',  and  z'  system  moves  through  space  with  the  payload,  it  can  be 
considered  tc  coincide  with  Johnston  It  land  if  only  rotation  between 
the  x’,  y',  z'  and  x,  y,  z  systems  is  to  be  measured  (Figure  2.1  (c) ). 
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Here,  A  and  ?  are  unit  vectors  pointing  in  the  direction  of  the  z'  axis 
(the  vehicle  longitudinal  axis)  and  the  earth's  magnetic  field  vector 

•f 

F,  respectively.  The  problem  la  to  determine  8  (t)  and  if  (t) ,  given 
8',  if',  A,  1  and  the  x',  y'  magnetometer  readings. 

The  polar  angles  6  and  f  are  given  by 


if 

e 


tan  ^ £  cos  (D,x)/  cos 
sin  1  ^  cos  (D,z)j 


(2.1) 

(2.2) 


where : 

cos  (D,x)  =  an  cos  e'  sin  ♦'  +  cos  0'  cos  if'  +  a31  sin  8' 

cos  (D.y)  =  a12  cos  8'  sin  if'  +  a22  cos  6’  cos  if'  +  aJ2  sin  0' 

cos  (D, z)  ■=  a13  cos  8'  sin  if'  +  a23  cos  0'  cos  if'  +  a33  sin  8' 

Here,  a^  are  the  transformation  (rotation)  coefficients,  i.e.,  direc¬ 
tion  cosines,  between  the  x  and  x'  systems.  Thus, 


x  *  anx  +  ai2y  +  ai3z 

y'  =  a21x  +  a22y  +  a2Jz  (2.3) 

z'  -  a31x  +  a32y  +  a33z 

The  vector  A  is  then  given  by 


A  =  i 


31 


+  j  a32  +  k  a 


33 


(2.4) 


-►  “►  ■+■ 

where  i,  j,  k  are  unit  vectors  along  the  positive  x,  y,  z  axes,  respec¬ 
tively.  Also, 


? 


i  cos  (F,x)  +  j  cos  (F,y)  +  k  cos  (F,z) 


(2.5) 


where  f  is  known.  The  determination  of  A  is  discussed  below.  The  re¬ 
maining  six  coefficients  ai2’  a13’  a21'  a22’  a23’  are  determined 

from  A,  ?  and  the  x',  y'  magnetometer  readings,  also  discussed  below. 
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2.2  DETERMINATION  OF  A 

2.2.1  General  Discussion.  Recall  that  A  is  a  unit  vector  along 
the  vehicle  longitudinal  axis,  i.e.,  the  vehicle  spin  axis.  As  the 
vehicle  accelerates  upward  from  the  launch  pad,  the  rush  of  air  past 
the  canted  fins  on  the  rear  causes  the  vehicle  to  spin.  As  the  vehicle 
gains  in  altitude  and  the  velocity  vector  V  begins  to  move  toward  the 
horizontal  (due  to  the  pull  of  gravity),  the  vehicle  will  tend  to  ac¬ 
quire  an  angle  of  attack,  i.e.,  the  velocity  vector  will  no  longer 
coincide  with  A,  but  will  lie  slightly  below  it.  If  the  vehicle  were 
not  spinning  (fins  not  canted)  it  would  feather  into  the  wind  like 
an  arrow,  i.e.,  A  would  rotate  with  the  velocity  vector,  pointing  in 
the  direction  of  V.  Since  the  vehicle  is  spinning,  the  torque  applied 
about  the  center  of  mass  by  the  air  drag  (the  drag  force  is  parallel 
and  opposite  to  V  and  is  applied  at  the  vehicle  center  of  pressure  , 
which  does  not  coincide  with  the  center  of  mass)  will  cause  the  vehicle 
to  precess  about  V,  which  is  referred  to  as  the  coning  motion.  If  the 
vehicle  is  low  in  the  atmosphere  and  if  V  is  rotating  rapidly  (a  low- 
summit  trajectory),  the  increasing  torque  will  cause  severe  coning  and 
will  eventually  cause  the  vehicle  to  spin  out  flat.  This  is  what 
happened  to  Rocket  18.  If  the  vehicle  is  climbing  rapidly  and  the 
direction  of  V  is  changing  slowly  (a  high-summit  trajectory),  the  cone 
angle  and  coning  frequency  will  remain  practically  constant,  and  the 
direction  of  the  axis  of  the  cone  will  remain  fixed  in  space,  coin¬ 
ciding  with  the  direction  of  V  at  the  instant  the  uniform  coning  motion 
is  established.  This  is  the  type  of  coning  eroerienced  by  Rockets 

15,  19,  and  26. 

If  the  initial  acceleration  is  large  and  the  launch  is  almost 
vertical,  the  drag  force,  and  hence  the  torque,  will  be  almost  non¬ 
existent  (due  to  low  air  density)  when  the  velocity  vector  finally 
begins  to  rotate,  and,  consequently,  the  vehicle  will  not  cone.  The 
direction  of  A  will  remain  fixed  in  space  coinciding  with  the  direc¬ 
tion  of  V  shortly  after  launch.  Rockets  8  and  9  did  not  experi¬ 
ence  any  perceptible  coning  motion. 
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Since  A  coincides  with  V  from  the  time  of  launch  to  the  time 
when  the  coning  motion  is  established,  or  for  a  short  time  after  launch 
in  the  case  of  Rockets  8  and  9,  it  is  important  to  know  V  during 
the  early  part  of  flight  in  order  to  determine  the  direction  of  A  for 
the  entire  flight.  The  computed  trajectories  (parabolic)  for  all 
rockets  begin  at  the  BRL  point,  which  occurs  at  varying  times  after 
launch.  (The  coning  motion  in  Rockets  15,  19,  and  26  begins  after 
the  BRL  point*  refer  to  Table  2.2.)  In  order  to  determine  V  from  the 
time  of  launch  to  the  BRL  point,  it  is  necessary  to  reconstruct  the 
trajectory  between  these  two  times.  This  can  done  by  a  cubic  fit, 
since  the  direction  of  V  is  known  at  launch  and  at  the  BRL  point.  Let 

■+  “f 

v  be  a  unit  vector  along  V,  i.e.. 


w 


-  i  v  +  j  V,  +  !c  v. 


(2.6) 


Then,  if  Bq,  i|>o>  and  0^,  are  the  azimuth  and  elevation  (refer  to 
Figure  2.2)  of  V  at  launch  and  the  BRL  point,  respectively,  the  direction 
cosines  are  given  as  a  function  of  z  (the  vertical  direction  in  the 
Johnston  Island  coordinate  system)  as  follows: 


v 


v 


v 


1 

2 

3 


S1  ^  +  sl^  +  3 
s 2  (1  +  s^2  +  s 

2  2 
(1  +  s j*  +  s2) 


Vl,! 


-1/2 


(2.7) 


where; 

s1  “  +  2b2  z  +  3b ^  z2 

2 

s2  “  Cj  +  2c2  z  +  3Cj  z 
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b  =  cot  ip  sin  8 
1  oo 

b„  =  3x.  -  (cot  ij<  sin  B.  +  2  cot  p  sin  B  )z. 

^  -L  i  1  o  o _ 1 

2 

Z1 

b,  =  -2x.  +  (cot  i/i.  sin  6,  +  cot  i p  sin  B  )z, 

■J  1  1  1 _ o  o  1 

3 

Z1 

c,  =  cot  ip  cot  8 

1  oo 

c2  "  3y^  -  (cot  cos  8^+2  cot  \j/Q  cos 

2 

Z1 

c,  «  -2y  +  (cot  p.  cos  B,  +  cot  ip  cos  8  )  z, 

■J  1  11  o  o  1 

3 

zi 

and  xx>  y1(  are  the  coordinates  of  the  BRL  point  in  the  x,  y,  z 
system. 

2_-2.2  Nonconing  Motion.  From  amplified  telemetry  records  of 
the  z-magnetometer  readings,  it  is  possible  to  determine  the  times  at 
which  coning  motion  begins,  for  Rockets  15,  19,  26,  and  the  times 
when  relative  motion  between  A  and  the  earth's  magnetic  field  F  ceases 
(i.e.,  when  A  acquires  a  fixed  direction  in  space),  for  Rockets  8 
and  9.  (Recall,  the  z-magnetometer  lies  along  A.)  Since  F  remains 
practically  constant  (in  magnitude  and  direction)  over  most  of  the 
flight,  a  nonvarying  z-magnetometer  reading  is  construed  to  mean  a 
nonvarying  A.  Table  2.2  shows  the  above  times  for  the  five  rockets. 
All  times  are  measured  from  launch.  The  last  column  shows  the  length 
of  time  that  A  coincides  with  v,  i.e.,  an  =  v^  a32  =  v2,  a33  =  v3 . 
During  the  period  from  launch  to  the  BRL  time,  v^,  v2 ,  v3  are  given  by 
Equation  2.7.  During  the  period  from  the  BRL  time  to  the  time  shown  in  the 
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last  column,  v^,  are  given  by  the  parabolic  trajectory.  For 

Rockets  8  and  9,  A  maintains  (until  reentry)  the  direction  it 
acquired  at  the  time  shown  in  the  last  column. 

2.2,3  Coning  Motion  of  A.  The  coning 

motion  of  Rockets  15,  19,  26  can  be  divided  into  two  phases:  coning 

buildup,  and  uniform  coning.  The  times  shown  in  the  third  column  of 

Table  2.2  actually  represent  the  beginning  of  the  coning  buildup  phase. 

During  this  phase,  the  amplitude  and  period  of  the  oscillatory 

z-magnetometer  (amplified)  readings  increase  from  zero  to  their  final 

constant  values  for  uniform  coning  motion.  In  order  to  determine  the 

vector  A,  we  then  proceed  as  follows.  Let  v  denote  the  unit  vector 

o 

which  lies  along  the  axis  of  the  circular  cone  (in  the  direction  of 
the  velocity  vector)  which  is  described  by  A  during  the  uniform  coning 
phase.  We  write 


-  iv 


10 


+  jv 


20 


+  kv 


30 


(2.8) 


The  method  of  computing  the  components  of  vq  is  described  later  in 

this  section.  We  now  construct  a  set  of  orthogonal  unit  vectors  e. , 
—♦  — #  1 

as  follows; 


f  -  <vo  •  f>VQ 

[i  •  <;„  ■  i)2),/2 


e3  X  C1 


(2.9) 


Refer  to  Figure  2.3.  Thus,  e1  lies  in  the  plane  defined  by  v~  and  f .  In 
Equation  2.9,  f  is  given  by  Equation  2.5,  and  70  by  Equation  2.8.  Now,  for 
circular  coning  motion 


A  =  sin  y  cos(wt-6c)e1  +  sin  y  sin(out-6c )e2  +  cos  y  e3  _ (2.10) 


where,  y  is  the  cone  half-angle,  U)  is  the  coning  frequency,  and  6  is 
the  coning  phase  angle.  Since  f  and  vq  are  expressed  in  terms  of  the 
Johnston  Island  coordinate  system,  therefore,  e^,  e2>  will  also  be 
expressed  in  this  system.  On  substituting  Equation  2.9  into  Equation  2.10,  A 
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will  be  given  In  this  system  and,  hence,  the  quantities  ^(t) >  a32^t^* 
a33(t)  can  be  readily  identified  (refer  to  Equation  2.4).  The  mechanics  of 
evaluating  e^,  e^,  e3 ,  substituting  them  in  Equation  2.10  and  identifying 
the  a3i(t) ,  can  be  programmed  on  the  digital  computer. 

In  order  to  determine  a33>  a^ ,  a33,  we  must  know  the  quantities 
Y,  w,  6^.  These  are  discussed  below. 

Uniform  Coning.  The  amplified  telemetry  records  of  the 
z-magnetometer  readings  for  Rockets  15,  19,  26  verify  that  the  vector  A 
ultimately  describes  circular,  coning  motion.  In  order  to  make  an 
accurate  determination  of  the  cone  half-angle,  y,  we  proceed  as  follows: 
Figure  2.4  shows  the  extreme  locations  of  the  coning  z'  axis  (z- 
magnetometer  axis)  which  are  coplanar  with  the  earth's  magnetic  field 
vector  F,  and  the  corresponding  locations  of  the  x'  axis  (x-magnetometer 
axis),  coplanar  with  F,  z',  z',  which  yield  the  maximum  positive 
x-magnetometer  readings,  Mxl  and  M  (recall,  x'  rotates  about  z'). 

Figure  2.5  shows  a  general  plot  of  the  x-magnetometer  reading  as  it 
appears  on  Rockets  15,  19,  26  during  the  uniform  coning  phase.  From 
Figure  2.4 

cos (tt/2  -  a.  )  =  M  ,/F  «  sin  o, 

1  xl  1 

cos(tt/2  -  a2)  =  Mx2/F  =  sin  <>2 

Hence , 

2y  =  or2  -  «1  =  sin-1(Mx2/F)  -  sin_1(M  j/F) 

where,  F  is  the  theoretical  field.  Although  the  above  equation  pro¬ 
vides  an  accurate  determination  of  y,  this  equation  cannot  be  employed, 
due  to  the  fact  that  y  is  too  small  (<  8°)  to  permit  an  accurate  and 

reliable  distinction  between  M  ,  and  M 

xl  x2 

The  relatively  small  coning  angles  experienced  on 
Flights  15,  19,  26  permit,  however,  a  simple  arid  accurate  computation 
of  y .From  Figure  2.6(a), cos  a  =  M  /F.  Differentiating  gives 

-sin  a  ■  tot  =  AM  /F 
z 
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where  £H^  is  defined  in  Figure  2.6(b),  which  shows  a  time  plot  of  the 

z-magnetometer  reading,  .  Since  the  coning  angle  is  small,  the 

differential  ior  can  be  written,  Ao  »  2y.  If  x',  Figure  2.6(a),  is 

coplanar  with  z  and  F,  then  the  x- magnetometer  will  read  a  maximum 

positive  value,  M  .  Hence, 
xmax  ’ 

cos(ir/2  -  a)  »  M  /F  ■  sin  a 
xmax 

Eliminating  sin  a  from  the  above  two  equations  and  disregarding  the 
minus  sign,  we  get 


Aa  -  2y  -  Z«z/Mxmax  (2.11) 

This  equation  has  the  advantage  of  not  containing  the  magnetic  field 
strength  F. 

The  coning  frequency,  u,  is  obtained  simply  by 


u  =  2  jt/T  (2.12) 

where,  T  »  period  of  uniform  coning. 

Since  f  (Figure  2.3)  lies  in  the  plane  of  vq  and  e^,  the 

angle  between  A  and  f  will  be  a  minimum  when  the  e^  component  of  A  is 

a  maximum.  However ,  when  this  angle  is  a  minimum,  the  z-magnetometer 

reading  will  be  a  maximum.  Let  t^  be  the  time  when  a  maximum 

z-magnetometer  reading  occurs  during  the  uniform  coning  phase.  Then 

the  e^  component  of  A  will  be  maximum  (Equation  2.10  )  when 

ut/  -  6  =  0,  or 

c 

6c  -  ut/  (2.13) 

In  passing,  we  point  out  that,  if  the  angular  momentum  vector  of 
coning  points  in  the  direction  of  v  ,  then  ui  >  0.  If  it  points  in  the 
opposite  direction,  then  u  <  0.  The  vehicle  cones  in  the  same  direc¬ 
tion  that  it  spins. 


Coning  Buildup.  If  we  replace  the  angle  lot  -  6c  with  the  symbol, 
X,  then  Equation  2.10  becomes 


A  -  sin  y  cos  x  ej  +  sin  y  sin  x  e2  +  cos  Y  e3 


(2.14) 
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where  x  is  defined  in  Figure  2.7.  During  the  coning  buildup  phase,  the 

angle  x  is  not  a  linear  function  of  tine,  due  to  the  fact  that  the 

coning  frequency  te  is  not  a  constant.  A  satisfactory  description  of 

X  during  the  buildup  phase  is  obtained  as  follows.  From  Figure  2.7  we 

see  that  the  z-magnetoraeter  reading  will  be  a  maximum  when  A  lies  in 

the  e, ,  v  plane  between  e,  and  v  ,  i.e.,  when  X  “  0,  and  that  it  will 
lo  1  o 

be  a  minimum  when  A  lies  in  this  same  plane  outside  of  e^  and  v  ,  i.e., 

when  X  =  tt.  Since  the  times  at  which  the  extreme  values  of  M  occur 

z 

are  known,  it  is  possible  to  make  a  graph  of  x  versus  t  according  to  the 
following  rule: 

X  =  2nn  (n  =  0,  1,  2....),  for  maximum  M 

z 

X  “  (2n+l)n  (n  •  0,  1,  2....),  for  minimum  M  (2.15) 

During  coning  buildup,  the  angle  y  increases  from  zero  to  its 
final,  constant  value  for  uniform  coning.  In  order  to  determine  the 
behavior  of  y  during  coning  buildup,  we  proceed  as  follows.  From 
Figure  2.7  we  see  that 

f  =  sin  X  e1  +  cos  X  e^  (2.16) 

Since  f  .  A  =  (WF,  where,  F  =  magnetic  field  strength,  therefore,  from 
Equations  2.13  and  2.16,  we  get 

Mz/F  =  sin  X  sin  y  cos  X  +  cos  X  cos  y 

Employing  Equation  2.15  we  then  get 

M  =  F(  sin  X  sin  y  +  cos  X  cos  Y) 
zmax 

M  .  *  F(-sin  \  sin  V  +  cos  \  cos  v) 

zmin 

From  the  above  two  expressions  we  readily  get 


whe  re , 


y  = 


tan 


-1 


cot  X  \ 


M  -  M  . 
zmax  zmin 

M 

zmax 


01  >  0) 


(2.17) 
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Before  Equation  2.17  can  be  employed  to  determine  y ,  we  must  first 
evaluate  X.  During  the  uniform  coning  phase,  we  select  any  consecutive 
values  of  M.^  and  Mzmin  and  evaluate  1).  We  then  calculate  the  value 
of  y  for  the  uniform  coning  phase,  which  is  obtained  from  Equation  2.11,  and 
substitute  y  and  1]  in  Equation  2.17  and  thereby  determine  X.  Thus, 


cot 


-1 


01- 


Note  that,  0  <  X  <  n . 


11 

Also,  M 


1  tan 


n  M 

TJ  zmax 


2  M 


xmax 


i] 


(2.18) 


zmax 


may  be  positive  or  negative,  but 


M  >0. 
xmax 

Having  evaluated  X,  we  can  now  determine  y(t)  from  Equation  2.17.  Thus, 

we  select  consecutive  values  M  ,  M  .  (or,  M  M  )  in  the 

zmax’  zmin  '■or’  zmin’  zmax-*  ln  the 
coning  buildup  phase,  evaluate  T|,  and  then  compute  y.  The  computed 

value  of  y  corresponds  to  the  instant  of  time  midway  between  M 
and  Mzmin ‘  By  all  the  and  readings  in  thT* 

buildup  phase,  we  can  therefore  determine  y(t)  throughout  this  phase. 

It  is  important  to  note  that  the  X  appearing  in  Equation  2.17  remains 
unchanged  during  the  computation  of  y(t),  i.e.,  it  is  the  single  value 
computed  from  Equation  2.18. 


Determination  of  v 


The  components  of  vq  are  the  final  quantities 


which  must  be  evaluated  before  A  can  be  determined  from  Equations  2.9  and  2.10. 
In  Figure  2.8  v0  and  f  are  defined  by  polar  coordinates  in  the  Johnston  Island 
coordinate  system.  From  Equation  2.8  and  Figure  2.8  we  see  that 


10 


sinu^  sindj 


v20  =  sinP’i 
v30  =  cos^l 


(2.19) 


The  angle  5j  is  nothing  more  than  the  azimuth  of  the  BRL  parabolic 
trajectory,  which  is  known.  We  must  then  determine  p  .  Since, 


f  •  v  =  cos  X 
o 

therefore  have 


v10  sintl2  sin52  +  v20  sinkl2  cos42  +  v30  cosP'2’  we 
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cos\ 


(2.20) 


sinn1  sinp2  cosfSg-S^  +  cosp^  cosp2  = 

The  A  appearing  in  Equation  2.20  is  that  computed  from  Equation  2.18.  At  John¬ 
ston  Island,  P2  ~  119.85“  ,  =  10.53“.  We  may  then  solve  Equation  2.20,  by 

graphical  means,  for  the  unknown,  p^. 

Coning  Characteristics  of  Rockets  15.  19.  26.  Table  2.3  shows  the 
salient  coning  characteristics  of  Rockets  15,  19,  26.  These  include 
the  duration  of  the  coning  buildup  phase,  the  angle  A  between  the  cone 
axis  vo  and  the  earth's  magnetic  field,  f,  the  polar  coordinates  (ji  6j) 
of  Vo,  the  cone  half-angle  y,  and  the  coning  period,  T(=  2tt/iju).  In  the 
case  of  Rockets  8  and  9,  vq  denotes  the  fixed  direction  in  space  assumed 
by  the  vehicle  longitudinal  axis,  A.  This  direction  is  simply  that  of 
the  vehicle  velocity  vector  at  110  seconds  after  launch.  Also  shown 
in  Table  2.3  is  the  approximate  frequency  of  spin  of  each  rocket  about 
its  longitudinal  axis,  which  was  determined  from  the  x-magnetometer  data. 

The  rather  unorthodox  value  of  p  computed  for  Rocket  19  is  due 
to  the  questionable  value  of  6^  employed  in  the  calculation.  Refer 
to  Section  3.2  for  a  discussion  <..  •  the  azimuth  of  the  Rocket  19  tra¬ 
jectory. 

Figures  2.9  through  2.11  show  y(t)  and  x(0  during  the  coning 
buildup  phase,  for  Rockets  15,  19,  26.  Note  that  the  slope  of  y(t) 
becomes  constant  toward  the  end  of  the  buildup  period,  indicating 
that  dx/dt  ■=(!)  =  constant  . 
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2.3  DETERMINATION  OF  ajj  ,  a12  »  a13  <  a21 1  a22  »  AND 

2.3.1  General  Discussion.  The  following  method  for  computing 
the  remaining  six  coefficients  a^,  a^,  ai$>  a2-,  a22’  a23  aPP^es 

— *  J 

equally  to  the  coning  buildup  and  uniform  coning  phases.  Let  i 

— *  f  f  f 

and  j  be  unit  vectors  along  the  x  ,  y  axes,  respectively,  i.e., 
along  the  sensitive  axes  of  the  x  and  y  magnetometers 
(Figure  2.1).  Thus 

i'  -  a11i  +  a12j  +  a1;}k 

(2.21) 

j  =  a21i  +  a22j  +  a23k 

Then,  the  unit  vectors  j',  A  form  an  orthogonal  set  (recall,  A  lieB 
along  the  z-magnetometer  axis).  We  now  construct  a  set  of  orthogonal 
unit  vectors  Uj^,  u2,  u3  as  follows  (Figure  2.12): 


U1  =  A>  u2 


.&:■&;  f]LA. 


[1-(A-?)2] 


1/2 


u3  =  u^  x  u2 


(2.22) 


During  coning  buildup,  A  Is  given  by  Equation  2.14,  and  during  uniform  coning, 
it  is  given  by  Equation  2.10.  From  Figure  2.4  we  see  that 


1^  =  cosa  u2  -  sina  u3 


j  =  sina  u2  +  cosa  u3 


(2.23) 


where; 


2  2 
(M  +  M  ; 
x  y 


1/2 


sina 


2  2 

(IT  +  M  ; 
x  y 


1/2 


(2.24) 
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Here,  M  and  M  are  the  x  and  y  magnetometer  readings,  respectively. 

— .  y 

Since  A  and  f  are  expressed  in  the  Johnston  Island  system,  therefore 
U'l' ,  u3  be  expressed  in  this  system  and,  consequently,  so  will 

i  and  J  .  By  comparing  Equations  2.23  and  2.21  it  should  then  be  possible 
to  identify  the  a ^ ,  a^,,  ai3>  a2l>  a22’  a23‘  mechanics  °f  eval¬ 

uating  u^,  u^,  u^,  i/,  i'  and  identifying  the  can  be  programmed  on 
the  digital  computer. 


2.3.2  Rocket  19.  Of  the  flights  being  considered  in  this  study, 
Rocket  19  is  unique  in  that  it  was  the  onlv  rocket  actually  in  flight  when 
burst  occurred.  One  consequence  of  the  burst  was  that  the  earth's  local 
magnetic  field  experienced  a  temporary  distortion,  which  lasted  for 
approximately  100  seconds.  During  this  time,  the  x,y,z  magnetometer 
readings  suffered  a  corresponding  perturbation,  thereby  making  them  un¬ 
suitable  for  the  purpose  of  attitude  determination.  It  became  necessary 
then  to  artificially  simulate  unperturbed  magnetometer  readings  M^  My, 
Mz  in  order  that  the  rocket  attitude  could  be  successfully  determined. 
Thus,  the  newly  generated  and  would  be  substituted  in  Equation  2.24, 
and  the  simulated  Mz  would  be  employed  to  ascertain  the  cone  half-angle, 
V ,  during  uniform  coning.  The  natural  Mz  could  be  employed  in  Equations  2.17 
and  2.18,  since  the  burst  occurred  after  the  start  of  the  uniform  coning 
phase . 

The  task  of  generating  M^,  My,  Mz  is  simplified  by  the  fact 
that  the  earth's  (unperturbed)  magnetic  field  changes  by  a  negligible 
amount  over  that  part  of  the  trajectory  in  which  che  perturbation  occurs. 
Thus,  the  earth's  field  can  be  regarded  as  a  constant.  The  method  of 
generating  MJC,  M  ,  Mz  is  as  follows.  In  Figure  2.13,  the  unit  vectors  p^, 

,  ^3  ^ornl  a  tight -handed ,  orthogonal  system  where 


=  cos(x  +  x/2)  +  sin(x  +  x/2) 

P2  -  Ax  p: 

P3  =  1 


(2.25) 
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Here,  x  =  wt  -  5c,  and  A  is  given  by  Equation  2.10,  where  to  =  constant. 
From  Figure  2.13  we  see  that 


i/  =  cosC  +  sinC  P2 
j7  =  -sin£  px  +  cos£  p2 


(2.26) 


Since  the  angles  x  +  */2,  y,  Q  are  the  usual  Euler  angles,  therefore, 
d£/dt  =  n  ■  d(x  +  it/2)/dt  +  cosy,  or  Q  m  Cl  -  to  cosy,  where  Cl  “  frequency 
of  spin  of  the  vehicle  about  the  A  axis.  Hence, 


C  =  (Cl  -  to  cosy)t  -  (2.27) 

where,  6 g  =  spin  phase  angle.  The  field  components  along  the  magnetometer 
axes  are  given  by 

Mx  =  FT(*  ‘  My  =  ’  )'),  Mz  “  FT(f  •  A)  (2.28) 

where,  F^,  is  the  magnitude  of  the  earth's  theoretical  magnetic  field, 

and  f  is  given  by  Equation  2.16.  Actually,  the  quantity  Fx  is  not  needed, 

since  it  cancels  out  when  M  and  M  are  substituted  in  Equation  2.24.  Also, 

x  y 

in  Equations  2.17  and  2.18,  Mz  occurs  as  a  ratio,  therefore  making  the  use  of 
Ft  unnecessary. 

In  order  to  determine  the  spin  phase  angle  6^,  a  value  first  must 
be  assigned  to  the  quantity  Cl  -  to  cosy,  which  is  the  rate  of  change  of 
the  angle  £  in  Figure  2.13.  Although  y  and  to  can  be  determined  from  the 
unperturbed  x  and  z  magnetometer  data  (recall,  the  burst  occurs  after 
the  start  of  uniform  coning)  ,  there  is  no  direct  way  to  evaluate  Cl,  the 
spin  frequency  of  the  vehicle  about  its  longitudinal  axis.  The  frequency 
of  the  oscillating  x  or  y  magnetometer  readings  is  not  equal  to  (1,  due 
to  the  presence  of  coning,  nor  is  it  equal  to  Cl  -  w  cosy,  since  the  fre¬ 
quency  Cl  -  to  cosy  represents  the  angular  rate  between  the  i'  direction 
(x  magnetometer)  and  the  moving  line  of -nodes,  p^  (Figure  2.13),  whereas 
the  frequency  of  the  x  magnetometer  represents  the  angular  rate  between 
i'  and  the  fixed  direction  7.  It  is  possible,  however,  to  evaluate 
Cl  -  to  cosy  by  making  use  of  Equation  2.28.  In  doing  so,  though,  we  assume 
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that  the  vehicle  indeed  describes  uniform  coning  motion  and  that  the 
actual  x,  y  and  z  magnetometer  readings  are  accurately  described  by 
Equation  2.28. 

Let  t;  be  the  time,  during  uniform  coning  (just  prior  to  burst)^ 
when  a  maximum  z  magnetometer  reading  occurs.  Then,  making  use  of  Equa¬ 
tion  2.13,  Mx  and  My  in  Equation  2.28  provide  us  with  the  following  relation: 

,  M 

6  »  (fl  -  w  cosy)t/  -  tan"1  (r~)  ,  (2.29) 

8  My  t 

Let  t' '  be  the  time  (during  uniform  coning)  when  the  x  magnetometer 

■>  ■+  / 

reading  passes  through  zero.  Then,  Equation  2.29  and  the  equation,  0  =  f  .  i  , 
combine  to  yield  the  relation 


fi  -  to  cosy  = 


— - —  f  tan 
it  _ i  L 


-1 


sinoif  t / ' -t ') 


t  -t 


cot\  siny  -  cosy  cosu(t,,  -tf) 


i  M  1 

-tan"1  (ji)  ,  t  (n-l)nl 


(2.30) 


Here,  the  integer  n  (=  1,2,3----)  denotes  the  number  of  zero  crossings 
of  M^,  following  t' ,  at  which  t**  is  taken.  Thus,  if  t' '  occurs  at  the 
first  zero  crossing  of  following  t  ,  then  u  =  1.  The  appropriate 
sign  is  chosen  so  as  to  make  Q  -  w  cosy  >  0. 

If  the  coning  motion  were  indeed  uniform,  and  the  spin  fre¬ 
quency  of  the  vehicle  about  its  longitudinal  axis  were  constant,  then 
we  would  expect  the  quantity  Q  -  u  cosy  to  be  a  constant,  independent 
of  the  choice  of  t'  or  t" .  Unfortunately,  however,  this  is  not  the 
case.  Table  2.4  shows  the  computed  values  of  fl-w  cosy,  from  Equation 
2.30,  for  five  values  of  t".  In  all  cases,  t*  =  107.6121  sec  (t  =  0  corre¬ 
sponds  to  launch) . 
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From  the  above  nonuniformity  of  Cl  -  U)  cosy,  we  must  conclude 
that,  either  the  coning  motion,  just  prior  to  burst,  is  not  actually 
uniform  and  that  fl  is  not  constant,  or  that  it  is  impossible  to  ac¬ 
curately  locate  the  times  t'  and  t "  from  the  digital  printouts  of 
the  x  and  y  magnetometer  data. 

It  is  possible  to  make  an  independent  evaluation  of  v 
(•  0  -  U)  cosy)  by  means  of  the  data  from  the  vertical  gamma  scanner 
detector.  If  both  the  vehicle  and  source  of  gamma  rayB  (burst)  were 
stationary  in  space,  then  the  frequency  f^  of  the  gamma  scanner  data 
would  equal  the  frequency  of  the  computed  azimuth  angle  cp  (Figure 
2.1  (a)  and  Equation  2.1)  of  the  gamma  scanner  (f  «  dcp/dt).  Since  the 
vehicle  is  in  motion,  the  quantity  f  -  f^  will  no  longer  be  zero, 
but  will  vary  from  point  to  point  along  the  trajectory.  It  is  pos¬ 
sible  to  precompute  the  value  of  f  -  fg  at  any  point  on  the  trajectory, 
since  the  azimuth  of  the  trajectory  and  the  coordinates  of  the  burst 

are  known.  Thus,  we  could  compute  f  -  f  at  some  specified  time,  t 

y  a  r  ’  ’ 

and  then  employ  that  value  of  \>  which  produces  the  appropriate  machine- 

computed  value  of  f  at  time  t,  i.e.,  the  f  which  makes  f  -  f  equal 
a  a  y  a  n 

to  the  precomputed  value.  We  would  then  regard  this  value  of  p  as 
the  correct  one. 

Unfortunately,  it  is  not  possible  to  employ  the  above  scheme  for 
determining  u,  since  the  azimuth  of  the  Rocket  19  trajectory  is  not 
well  known  (Reference  1,  page  232). 

As  a  consequence  of 

this  uncertainty  in  the  trajectory,  it  was  decided  not  to  conduct  any 
data  reduction  on  Rocket  19  at  this  time.  Presumably,  if  a  more 
reliable  estimate  of  the  trajectory  azimuth  could  be  obtained,  it 
would  be  possible  to  complete  the  analysis  of  this  particular  flight. 
This  problem  is  presently  being  pursued  under  the  3850  data  analysis 
program, since  the  determination  of  the  parameters  which  may  indicate 
the  proper  trajectory  lies  in  an  analysis  of  the  Rocket  19  data. 
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2.4  ROCKET  18 

Recall  that  the  ultimate  purpose  of  this  study  is  to  describe  the 
shape,  size,  and  nature  of  the  debris  cloud  on  the  basis  of  data  taken 
from  the  photometer,  p  -detector,  and  the  four  gamma  scanners.  This  will 
be  done  by  comparing  the  observed  data  from  these  six  detectors  with 
the  computed  (as  a  function  of  time)  attitude  of  the  detector  axes. 

This  task  will  be  relatively  simple  in  the  case  of  Rockets  8,  9,  15, 
and  26,  since  these  payloads  are  spinning  at  a  reasonably  constant  rate. 
However,  in  the  case  of  Rocket  18,  the  task  would  be  formidable,  since 
the  spin  rate  of  this  rocket  is  extremely  erratic,  as  evidenced  by  the 
x  and  y  magnetometer  readings.  According  to  the  telemetry  records,  the 
payload  does  not  exhibit  a  well-defined  spin,  but  instead  appears  to 
rotate  aimlessly  about  the  longitudinal  axis  (the  longitudinal  axis, 
however,  describes  a  uniform  90*degree  coning  motion).  Due  to  the 
irregular  behavior  of  the  spin  frequency,  it  was  decided  not  to  perform 
a  data  reduction  study  on  Rocket  18. 

2.5  ROCKET  29 

Attitude  determination  was  not  performed  on  this  rocket,  due  to  the 
fact  that  data  were  low  level  and  unreliable. 
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TABLE  2.1  POINTING  DIRECTION  OF  DETECTOR  AXES 


Detector 

e’ 

♦  ' 

Photometer 

0 

906 

8-detector 

0 

90° 

Horizontal 

gamma  scanner 

-15° 

45° 

Vertical 

gamma  scanner 

0 

0 

TABLE  2.2  TIME  PARAMETERS  FOR  CONING 


Rocket 

BRL  Time 

Beginning 

of 

Coning 

When  A  assumes  a 
fixed  airection 

A  =  v,  from  launch 
until  -  -  - 

sec 

sec 

sec 

sec 

8 

110 

— 

110 

110 

9 

110 

— 

11C 

110 

15 

40 

45.6 

— 

45.6 

19 

30 

48.6 

.... 

48.6 

26 

35 

31.0 

... 

31.0 
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TABLE  2.3  CONING  CHARACTERISTICS 


Rocket 

Duration 
of  Coning 
Buildup 

\ 

V 

0 

Uniform  Coning 

Spin 

Frequency 

(cps) 

“l 

61 

_ X _ 

T 

8 

— 

— 

30°27 1 

26°12  ’ 

1 

... 

6.10 

9 

— 

— 

21°39 ' 

23°30 ' 

... 

... 

6.55 

15 

69  .37  sec 

93°40' 

26°36 1 

21° 

6°35 ' 

25.75  sec 

1  .52 

19 

52.29  sec 

113°52 ' 

-10° 

135°(?  ) 

7°56  ' 

25.70  sec 

2.05 

26 

77.75  sec 

121°  4' 

6°17 ' 

113°30' 

1°31 ' 

22.00  sec 

2.44 

TABLE  2.4  n  -  w  COSy  FOR  SEVERAL  VALUES  OF  t " 


t"-t' 

0  -  p)  cosy 

_ 

0.1076 

sec 

12.42611 

rad/ sec 

+ 

0.1379 

sec 

12.67248 

rad/sec 

+ 

0.3809 

sec 

12.68755 

rad/sec 

+ 

0.6259 

sec 

12.64988 

rad/sec 

+4- 

5.7526 

sec 

12.56068 

rad/ sec 
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Z  (Vertical) 


(0) 


2  IWagmrjd'ile'! 


(Magnetometer) 


Z 


Figure  2.1  Definition  of  Johnston  Island  and  payload  coordinate  systems. 
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Figure  2.6  Derivation  of  y  and  definition  of  AM; 


Figure  2.7  Definition  of  \  and  x  • 


Z  (Vertical) 


Y  (North) 


X(Ea«t) 


Figure  2.8  Definition  of  m ,  M2  >  6i  .  and  $2 . 
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et  26. 
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CHAPTER  3 

COMPUTING  PROGRAM  DESCRIPTION 


3 . 1  INTRODUCTION 

The  attitude  determination  program  described  in  this  chapter  was 
designed  through  modification  of  a  similar  program  constructed  by 
General  Electric  for  the  same  purpose.  The  program  computes  detector 
directions  in  the  manner  developed  in  the  theoretical  discussion  of 
attitude  determination  in  the  previous  chapter,  and  the  mathematical 
symbols  and  computer  symbols  are  related  wherever  possible. 

The  data  tapes  for  this  program  were  produced  by  General  Electric 
by  analog-to-digital  conversion  of  telemetry  data.  Considerable 
patching  of  the  tapes  was  apparently  necessary  to  eliminate  noise 
from  the  digital  information.  The  discontinuities  produced  by  this 
patching  required  special  data-smoothing  techniques  to  avoid  smooth¬ 
ing  continuously  over  the  discontinuities.  This  smoothing  process  was 
performed  on  all  digital  tapes,  and  computer  program  data  input  tapes 
were  produced  for  each  rocket. 

The  data  processing  with  the  attitude  determination  program  was 
done  at  Service  Bureau  Corporation  in  El  Segundo,  California,  on  an 
IBM  7094  computer.  Three  sets  of  data  printout  and  one  input  tape 
for  the  3850  Data  Analysis  Program  were  produced  for  Rockets  8,  9,  15, 
and  26.  The  data  from  the  remaining  Rockets  18,  19,  and  29  were  not 
processed  for  the  reasons  outlined  in  the  previous  chapter. 

The  processed  data  are  unclassified;  however,  certain  classified 
inputs  were  necessary  to  the  computer  program.  These  inputs  are 
discussed  in  Reference  2. 
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3.2  DATA  SMOOTHING  PROGRAM 

The  smoothing  program  consists  of  a  main  control  program  with 
associated  subprograms  that  will  smooth  data  on  multiple  tapes.  The 
program  checks  for  time  discontinuities  and  makes  appropriate  changes 
in  the  procedure  for  removing  bias  and  in  the  use  of  the  appropriate 
smoothing  formulas. 

Separate  routines  were  written  for  input,  output,  data  scanning, 
actual  smoothing,  and  bias  determination. 

The  procedure  for  the  smoothing  of  the  data  was  a  third-degree 
seven-point  least  squares  method.  The  formulas  used  were  obtained 
from  Reference  3.  The  coefficients  of  the  seven  smoothing 
formulas  were  arranged  in  the  seven-  by  seven-square  matrix  A  (by 
subroutine  ARRAY),  and  then  the  smoothing  was  performed  by  forming  a 
vector  x  of  seven  functional  values  to  be  smoothed  and  multiplying 
this  vector  by  the  matrix  A  to  obtain  a  vector  y  of  smoothed  points, 
i.e., 

y  »  Ax 

The  center  smoothing  formula  was  used  on  all  points  that  were  at 
least  three  points  away  from  a  data  discontinuity.  After  the  smooth¬ 
ing  of  a  block  of  continuous  points,  the  bias  for  the  particular  block 
was  determined  and  removed,  and  the  block  was  then  output. 

3.2.1  Main  Control  Program.  This  program  is  written  in  Fortran  II 
for  the  IBM  7094  and  acts  as  a  monitor  for  several  subprograms.  The 
program  inputs  start  and  stop  times  for  smoothing  and  printing,  run 
numbers,  and  tape  reel  numbers.  In  order  to  accomplish  a  multiple- 
tape  input  capability  with  a  minimum  of  operator  setup  error,  on-line 
messages  are  given  to  the  machine  operator  concerning  where  and  when 
specific  tapes  are  to  be  mounted  for  input. 

3.2.2  Subroutine  SMOOTH.  This  routine  written  in  Fortran  II 
performs  the  actual  smoothing  of  the  data  and  transmits  the  smooth 
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data  to  the  output  subroutine.  The  routine  calls  subroutine  ARRAY  to 
initialize  the  smoothing  operator  and  then  uses  this  matrix  to  do  the 
smoothing.  The  matrix  A  consists  of  the  coefficients  in  the  following 
equations : 

yi-3  =  h  (39xi-3  +  8xi-2  -  4xi-l  -  4xi  +  xi+l  +  4xi+2  -  2xi+3> 

yi-2  =  h  (  8xi-3  +19xi-2  +16xi-l  +  6xi  ’  4xi+r  7xi+2  +  4xi+3> 

yi-l  =  h  ('4xi-3  +16xi-2  +19xi-l  +12xi  +  2xi+r  4xi+2  +  xi+3> 

yi  “IT  ('2xi-3  +  3xi-2  +  6xi-l  +  7xi  +  6xi+l+  3xi+2  *  2xiJ-3^ 

yi+l>  yi+2>  yi+3  . 

where  is  the  smooth  value  corresponding  to  x^,  which  is  the  center 
point  of  the  seven  points  considered  by  the  smoothing  values. 

The  first  three  equations  were  used  on  the  first  three  data  points 
following  a  discontinuity,  the  last  three  on  the  last  three  data  points 
before  a  discontinuity,  and  the  center  one  on  all  other  points. 

3.2.3  Subroutine  ARRAY.  This  routine  is  written  in  FAP  and  is 
executed  only  once  at  the  beginning  of  the  program  execution.  It 
initializes  the  smoothing  matrix  A  which  consists  of  the  coefficients 
of  the  smoothing  formulas. 

3.2.4  Subroutine  SCAN.  This  subroutine  written  in  Fortran  II 
calls  the  input  subroutine  and  scans  the  data  to  find  a  sequence  of 
continuous  points.  Once  it  finds  at  least  seven  consecutive  points 
it  calls  subroutine  SMOOTH  to  perform  the  actual  smoothing,  and  then 
scans  ahead  to  find  the  next  discontinuity. 

3.2.5  Subroutine  GET.  This  buffered  input  routine  is  written 
in  FAP  for  the  7094.  It  reads  in  a  block  of  200  data  vectors  into 

core  storage  and  makes  this  matrix  available  to  both  the  LOAD  subroutine 
and  the  SCAN  subroutine.  The  routine  also  checks  for  tape  reading 
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errors  and  makes  five  attempts  to  read  the  tape.  If  after  five  tries 
the  error  remains,  the  record  is  skipped  and  input  continues. 

-3-'2-6  .  Subroutine  OUTPUT.  This  Fortran  II  subroutine  prepares 
a  block  of  100  output  vectors  for  storing  on  binary  output  tape  and 
printing  on  BCD  print  tape.  It  determines  the  amount  of  bias  in  the 
smoothed  data  by  finding  the  minimum  and  maximum  points  on  the  x  and 
y  magnetometer  values  and  calling  subroutine  PARAB  to  pass  a  least 
squares  parabola  through  five  points  including  the  minimum  or  maximum 
point  as  the  center  one  and  finding  the  minimum  or  maximum  of  this 
curve.  Once  it  has  obtained  the  best  estimate  of  the  minimum  and 
maximum,  the  bias  is  then  removed  from  the  whole  output  block. 

— — ' ?  — Subroutine  ERROR.  This  subroutine  writes  BCD  error  message 

on  output  tape  6,  when  errors  occur  on  reading  the  binary  input  tapes. 

— Subroutine  PARAB.  This  subroutine  is  written  in  Fortran  II 

and  determines  the  minimum  or  maximum  of  a  least  squares  parabola  fitted 
to  five  points  which  contain  the  argument  point  as  its  center. 

— Subroutine  LOAD.  Subroutine  LOAD  is  written  in  Fortran  II. 
This  routine  covers  the  special  case  when  a  discontinuity  appears  within 
seven  points  of  the  last  point  in  an  input  array.  To  handle  this  situa¬ 
tion  subroutine  LOAD  will  call  subroutine  GET  to  input  a  new  block  of 
points  and  place  the  remaining  points  of  the  last  block  in  front  of  the 
new  block  and  return  control  to  subroutine  SCAN. 

~---°  ?lw  Diagrams,  The  following  diagrams  outline  the  logical 
and  computational  processes  for  each  of  the  routines  discussed  above 
and  are  followed  by  the  program  coding  as  stated  in  Fortran, 
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Writ*  F..O.F. 
On  Tip* 
Nod 


Writ*  On-Line 
Messog*  ToOporolor 
"Prtpor*  Top*  No  X 
A*  Th»  Next  Input 
Top*  On  B-5" 


Rewind  And 
Unload  Tap* 
No.  9  i 


tmr  SWrtiT. 
I >w3n 


Wtit*  On-LIn* 
Message  To  Operator 
Pr  pare  N*w  Output 
Save  Top«ForN*xt 
Output  Tap*  On  A-s" 


Writ*  On -Un* 
Message  To  Operator 
It  New  Output  Tap* 
Hoe  Been  Loaded 
On  A-5,  And  If 
Tap*  No.  X  it 
Ready,  Push  Slort* 


(CALL  SCAN) 


Wi  >■  Hun  Nm  S 
Or  Tjt*  M,  J 


R«inltiolii* 

I/P  Parameters 
IRun=II,ITPNO  =  I2 
STRTSM=X3,FNSHM  =X4 
STRTPR=  X5,FNSHPR=X6 
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46 


48 


49 


Calculate 

Minimum 
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_ £ _ WAIN  PROGRAM  FOR  DATA  SMOOTHING 

ID ( 24B4 ) BUFFER(  12'  2°0.  2)  .DATA  112,207), OUT  <15, 4661,8(2466,2), - 

. . . 

inn  gSggP^f  ^ilOO.tRUN.IIkNU.STBTSH.FUSHSW.STRTPft.FWiHPft - 

WRITE  TAPE  9,  IRUN  ' — - — — - — _ __ 

- 2-QQ-EM1NT  300,  ITPNO _ 

35?  ' 1H111HIF  TAPE  NO*I5'29”H  IS  REAOV-ON  B-5.  PUSH  SUrTH - 

PAUSE  * -  - - - - - - - - - - 

- SfAL(N^UJ_LAP^^lOO_LIAi.Ut.X3J  X4 f  X5 ,  X6 

IF(Il)  400,1100,400  - — - - - - 

— AM_£mT.  500.1?  _ _ 

500  IAPE  N°:r5’31‘H  A'S  T HE-nWInput— APE  ON  B-sT7~ 


600  I SW  1=2 


^PARE  NEW  OufpUnivrfAPE  FOR  THE  NEXT 


800  CALI  SCAN 


I TPNO=I 2 
STRTSH=X3 
FNSHSM=X4 
-5.IRTPR=X5. 
FNSHPR=X6 


900  ENO  FILE  9 


PRINT  1000, ITPNO  '  ~ — - - - 

'  «°r m  push^tart?!0^  °"  m 

WRITE  TAPE  9,  IRUN  ~  - - - - - - 

- SO  TO  35  0 _ _ 

1100  PRINT  1200  * - - - - - - - 

1 2-QQ  -P3PMAJ15 3H  IN 0_MQR E  _XA-P  £.S  _WXUL.  -B E_  NE E D  ED_, — IMIS_  J^_T^H^FJJ^J1_^UNj 
_ END  FILE  9 

CALL  REWULD  (9)  ’  — “ - — - - - 

_ CALL  E X  I  T _ 

END  "  .  — - - - - - - 
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_ 5UWM1UTJNE  SMOOTH  (  i  ,j,K ) 

MNJHSJON  BUFFER  ( 12*200. 2  I  .DAT  A  <12.207 )•  OUT  (15*200 )»B<24od*2). 

COMMON  BUFFER, DATA, OUT,  IRUN.STRTSM.FNSHSM.STRTPR.FNSHPR 
EQUIVALENCE  ( BUFFER ,B ),( DATA ,D ) 

DIMENSION  A (7,7) 

_  IFllllll  >10,20,10 

10  CALL  ARRAY(  11111  ,A)  . ' 

- 20  IF( K-4 ) 200, 300  *300 

200  MSTART-I-1 
_  GO  TO  400 

300  MST  ART  *  I  —4  - - 

400  DO  1000  L*7,9 

500  OUT ( L , J ) >0,  - - - 

600  DO  900  M«1 ,7 

700  MM«MSTART+M  . 

_  *00  OUT(L*J)*OUT(L.J)-»A(K,MI#DATA(L*MM> 

900  CONTINUE 
1000  CONTI N U E 

1100  I F ( K-7 11150,1300 ,1300  ‘  '  - - - 

_ 1150  IF(K-1 >1200,1300,1200 

B1200  OUT  (  1 5 ,  J )  *606060606060  . . . 

GO  TO  1400 

B 1 300  OUT ( 1 5 ,J ) “546060606060 
- L400  GO  10(1401 , 1401 ,1401 j 1402 , 1403,1403,1403 ) ,K 

1401  CALL  OUTPUT! t-l+K , J.K ,0 )  - - 

_  GO  TO  1500 

1402  CALL  OUTPUT (  I«J,K«0) 

GO  TO  1500 

1403  CALL  OUTPUT (I-4+K.J.K, 01 

_ 1500  RETURN _  _ 

END  .  '  - - - 
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PAP 

COUNT 

ENTRY 

ARRAY  STZ* 
CIA 
ADO 
STA 
SX  A 
AXT 

82  CLA 

82A  STO 
T I  X 

B3  A_XT 

TRA 

DATA  DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
DEC 
FND 


20 

ARRAY 

1.4 

2.4 
-1 
82A 
88,1 

49.1 

DATA+49,1 

**.l 

82.1.1 

**,1 

8.4 


.92857143,. 

.09523810,- 

.14285714,- 

.38095236,. 

.02380952,- 

.28571429, . 

.04761905,.; 

.09523810,- 

.45238095,. 

-.09523810,- 


1 19047619, - 
-.04761905, 
-.09523810, 

45238095.. 
-.09523810, 
142B57] 4  ,- 

28571429..  < 
.16666667,- 
19047619,-, 
-.09523810, 


.09523810, -.09523810,. 02380952 
.1904761 9,. 45238095,. 38095238 
-. 16666667, .09523 8 10, -.0952 38 10 

28571429. . 04761 905, -.09 52 38 10 

. 1428 571 4,. 28571429,, 333133J3 _ 

. 09523810,. 02 380952, -.09523810 

45238095..  38095238, -.09523810 
-. 095 2 3R 10. .14285714, .38095238 
.0476 1905,. 095238 10.. 02380952 

,.19047619, .92857143 


J 

f 


5 


I 
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SUBROUTINE  SCAN  ....  . . . 

DIMENSION  SUFFER (1 2, 200.2  It OAT A  1 12.207 ), OUT <15. 200  I  •  B ( 2400.2 ). 
10(2484) 

COMMON  BUFFER. DATA,OUT« IRUN.STRTSM.FNSHSM*STRTPR,ENSHPR 
EQUIVALENCE  ( BUFFER , B ) . ( DATA .D I 
100  CALL  GET! IB, EOF, ISIZE) 

- 110  IF[STRTSM-BUFFER(1, 200, 18) 1120*120*100  _ _  _ 

120  I F ( FNSH5M-BUFFER (1,1, 18)1130*130,140 
130  RETURN 
140  CONTINUE 

DO  150  1=1, ISIZE  ...  _ 

150  D( I)»B(  I, IB) 

_ J*1  .......  _  _ 

ISIZE*  ISIZF/12 
200  I«1 

300  IF( ISIZE-I-7)400,700,700 
400  IFIEOE >600,500,600 
500  CALL  OUTPUT  (I.J,K,1.) 

_ RETURN  _  _  _ 

600  CALL  LOAD (ISIZE, EOF) 

GO  TO  200 
700  ISTART=I 
I  END* I +5 

DO  740  I X* I  START , I  END 

_ TEST=DATA(1,IX+1)-DATA(1,IX>  _  _ _ _ 

710  IFITEST-. 005)750, 720, 720 

_  7.20  IFITEST-. 0151740, 740, 750  _ 

740  CONTINUE 
GO  TO  800 
750  I = I X+ 1 

_ GQ._TO  .300  . .  . . . 

800  DO  1000  K=1 , 3 
1000  CALL  S MOOT H  ( I « J ,  K  ) 

1100  1=1+2 

__1 2  0  0  1*1  +  1  _  _ 

1300  CALL  SMOOTH! I, J, 4) 

2. lACHL  I  Ft.  u;_i  Z  Erl -4)2.000  *  1500, 1500  .  ...  _ 

1500  TEST=DATA( 1,1+4  I -DATA! 1,1+3  ) 
tFJ.TE.ST-,  005  11600, 151 0,1510 
1510  IF( TEST-. 0] 5 11200,1 600, 1600 
1600  DO  1800  K*5,7 
1800  CALL  S  MOO  T  H  ( I ,  J « K ) 

_ 1900  I  *  1+4 

GO  TO  300 

2000  IFIEOF >2400,2100,2400 
2100  DO  2300  K=5 ,7 
2300  CALL  SMOOTH!  I  ,J,K) 

GO  TO  500 

_ 2400  CALL  LOAD! ISIZE, EOF)  _ 

1*4 

GO  TO  1400 
END 
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•GET  SUBROUTINE  FOR  INPUT 

_ COUNT  _  50 _ 

ENTRY  GET 
_ ENTRY 


ZERO  STZ*  1,4 

- ISA.  2.4 

GET  LMTM 


STZ  FIRST 

_ _ 

SXO  BUFNO , 5 

_ BXDfl _ i _ 

RCH8  SKIP 


•  +  1 
5 


INPUT, 5 
* 


CLA  BUFER2+SIZE6.6 

STO  BUFER2+SI Z El 2 , 7 

-SIS _ BUFER2+  S 1 2  Fft  ■  4 

TRA  B6C 


BUFNO,  5 
_ BUFER , 5 


STZ* 

2,4 

RUNS 

5 

STL 

FIRST 

LDC 

BUFNO,  5 

BUFNO ,  5 

IRA 

Bl  1 A 

B15 

SXA 

♦•*•2,4  . _ 

TSX 

$ERROR,4 

AXT _ 

♦  ♦,4 

TRA 

B4 

m  ■mini  mi  hi 

SIZE 

EQU 

,’na 

SIZE12 

_EQU_ 

FIRST 

PZE 

i 

SKIP _ IQCjl 

U.fl  .12 

IOCD 

BUFE  R 1 , 0  ,  S I Z  E 1 2 

IOCD 

BUFER2 ,0,SIZ£12 

BUFNO 

..EQU 

INPUT 

PZE 

0,0,1 

BUFFR 

PZE _ _ 

PZE 

0,0,2 

--SIZES 

__EjQU 

SIZE*6  _ 

TWO 

PZE 

0 , 0 , SI Z  E 12 

_ COUNT  -EQU _ 

FIRST  .  _ . 

BUFER1 

EQU 

32562-SIZE12*2 

_ BUFF&Z 

-EQU _ 

3?SK?-Sr7F12 

END 

WRITE  OUTPUT  TAPE  6,  1 

RETURN _ 

1 

FORMAT (5X27HERROR  IN  READING  INPUT  TAPE) 

■ 

END 
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-  SUBROUTINE  OUTPUT  ( I  l  »J.,ICK  .FLAG) 

DIMENSION  BUFFER (12. 200. 2  )  .DATA 1 1 2 .207 > .OUT  1 15 • 200 ) *B (2400.2 ) • 
10(2484) 

COMMON  BUFFER, DATA, OUT. I  RUN .STRTSM. FNSHSM * STRTPR .FNSHPR 
EQUIVALENCE  (BUFFER , B ),( DATA ,D I 
DIMENSION  MI  NMAX ( 4 ) , AMPL ( 4  ) 

-  EQUIVALENCE  (MINMAXI  1 )  •  1 XM I N ) ,  ( M INMAX (  2  )  ,  1 XMAX  )  ,  (MINMAXOl  .lY-MLNU 

1 ( M I NMAX ( 4 ) , I YMAX ) , ( AMPL ( 1 )  ,XM! N ) , ( AMPL ( 2 1 » XMAX ) , ( AMPL ( 3 ) , YM IN ) , 

2 ( AMPL (41# YMAX ) 

!eii 

IEND'100 

IFIlIUt  )  50,100,50  . 

_ 5Q..  X_GA  IN«1. 

YGAIN'1. 

BIAS'O. 

CALL  7ER0  (11111 ) 

100  IF(FLAG) 20 0,300, 200 

200  IFND'J  - 

_ 50  TO  600  _ 

300  DO  320  L=1 ,6  - 

320  OUT(L»J)'DATA(L«I) 

DO  350  L'10,12  - 

350  OUT (L.J)'DATA(L.I) 

400  J'J+1 

500  \F( j- 1  END  I  2700 t2700, 600 

600  IXMAX'l  - - 

_ I XM IN' 1 

IYMAX'l  •  . . — 

1 YM I N  = 1 

XMAX'OUT (7.1)  '  - - - 

_ .XMIN'OUTt  7,1) 

YMAX 'OUT (8,1)  — - 

YMIN'OUT (8,1) 

1 1600' IENO-1  ’  . . - 

700  po  1600  1=2,11600 

8^0  IF(OUT(7,I)-XMAX)1000, 1000. 900 
_ 9_0_0  J  XMAJt'I 

XMAX'OUT (7, 1)  . . 

GO  TO  1200 

1000  !F(OUT(7,I)-XM!N)1100, 1200, 1200  - 

1.100  IXMIN'I 

XMIN'OUT (7,1 )  -  - 

?00  ..I.FCOUT  («  » I )  -YMAX  )  1400 , 1400,1300 

1.700  IYMAX'l  - - 

YMAX'OUT (8,1) 

GO  TO  1600 

1400  I F( OUT ( 8, 1 l-YMIN ) 1500, 1600.1600 
1500  I YM IN= I 
_  _ YMIN'OUT (8, I) 

1600  continue  . . . 

1.7.00  DO  2000  K'1,4 

1800  CALL  PABA8 (MJNMAX , AMPL, K.IEND, TAG) 

IF(T AG  1220 0,20 00, 2 200 
2000  CONTINUE 
2100  XB I AS' • 5* ( XMAX+XMI N ) 

XMAX'XMAX-XBIAS 

YBIAS'.5»(YMAX+YMIN) 

YMAX'YMAX-YBIAS 

GA I N'SQRTF ( XMAX*YMAX ) 
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_  XGAIN'GAIN/XMAX. 

YGAIN=GAIN/YMAX 
RATTO*XMAX/YMAX 
2200  00  2400  1*1 , IEND 

OUT (7. 1)* (OUT (7. I ) -XB I  AS ) »XGA I N 
OUT  (8,1)  =  (OUT (8 . 1  >  — YB I  AS  )*YGAIN 

_ OUT (14. 1 ) '1 .- ( OUT (8.1 )  /YMAX I *»2- ( OUT  17.1 )/XMAX)**2 

2400  OUT( 13.1 )*RATIO 
C  OUTPUT  ROUTINE 

IF(STRTPR-OUT( 1 ,IENO)  12450.3050.3050 
2450  I F ( FNSHPR-OUT (1.1 1)3050,2500.2500 
2500  WRITE  OUTPUT  TAPE  6.2600. IRUN 

-?M0  FORMAT  (43H10UTPUT  FROM  DATA  SMOOTHING  PROGRAM  RUN.  NO,I.3»18HEOi  “ 
1M  RUB INSTE1N///6X4HTI ME 13X1 HY 15X1 HX 15x1 HZ  1 3X5HRATI 012X4HD1 FF8X4HCO 
2DF 1 

IF (50- I FND) 3100,2800,2800 
2800  DO  2ROO  KK*1 , IEND 

2900  WRITE  OUTPUT  TAPE  6* 3000 .OUT ( 1 »KK  > .OUT 1 7.KK )  .OUT ( 8.KK ) .OUT ( 9. KK )  .0 

_ 1UTJ  1 3 » KK  )  .  OUT  (  14,  KK)  »0UT(15.KK  ) 

3000  FORMAT(F12.4,F14.4,2F16.4,2F16.3»10XA1 ) 

3050  WRITE  TAPE  9 , 1  END , ( (OUT  1 1 , J ) . I *1 , 1 5  I  ,  J«1 ,  I  END ) 

J  =  1 

2700  RETURN 

3100  DO  3200  KK=1 .50 

- 3.2P0  WRITE  OUTPUT  TAPE  6 . 3000  .OUT  ( 1  «KK  )  .OUT  I7.KK  )  .OUKB.KK  ).0UT(9.XK).n 

1UT( 13.KK ) .OUT ( 14.KK ) »OUT( 15.KK) 

WRITE  OUTPUT  TAPE  6.2600, IRUN 
DO  3300  KK*51 , 1  END 

3300  WRITE  OUTPUT  TAPE  6 . 3000, OUT ( 1 .KK > .OUT ( 7.KK ) .OUT ( 8 .KK I .OUT  (  9 .KK I .0 
1UT (13, KK ) .OUT ( 14.KK) .OUT ( 1 5.KK ) 

. .  WRITE  TAPE  9,IEND,((0UT(I,J),I  =  1.15),J.1,IEND) 

1 

RETURN 

END 
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_ SUBROUTINE  PARAB  I.M1NMAX*  AMPL.K*  I  END*  TAG)  .  . . . 

DIMENSION  PUFFFRI 1 2, 200 . 2 > .DAT A ( 1 2 , 207 > .OUT f 1 5 *200 ) *B 1 2400. 2 ) . 

I D ( 2484  > 

COMMON  BUFFEP, DATA. OUT. IRUN 

EQUIVALENCE  (BUFFER, B) .( DATA, D»  . - 

DIMENSION  MINMAX ( 4 ) , AMPL (  4  I 

_ 100  I_»M  INMAX  (K)  .  .  ..  _ _ 

IF<  1-1)300,300,400 

_ 300  TAG=1.  . . . . 

GO  TO  1400 

400  IF(  I-IENP'500.300,300  . . . . 

R  500  TEST*6 06 06 0606060 

_ IF_(_0_UT  (15iL)-TFS_T  1300,600 ,300  _ 

600  TAG-0. 

700  I F (OUT (15 , 1+1 )-TEST) 000,900,000  _  ... . . . . 

800  1=1-1 

GO  TO  1300  _ _ 

900  1F(  I  +  1-IENDI1000, 800, 1000 

1000  I F(0UT(_15,I-lJ-TESmi00, 1200,1100.. _ _ 

1100  1=1+1 

GO  TO  1300  _  _ 

1200  IF(  1-211300,1100,1300 

1300  L=(13+K)/2  _ 

AO* ( -6 •* (OUT ( L , I -2 ) +OUT ( L  ,  1+2 ) )+24.* (OUT ( L , 1-1 l+OUT ( L » I +1 ) ) 

I  +  34.+0UT  ( L .  I  M/70.  _  .  ...  _ 

A1  =  ( 1 4 ,*(OUT  (L  » I  +2  )  —OUT  ( L  ,  1-2 )  >+7.*(0UT(l.  ,  1+1  )-OUT(L,  1-11)1/70, 
A2=(lO.*(OUT(L,I-2)-OUT(L.n+OUT(L,I+2)  ) -5.*( OUT (L,  1-1  J+OULILjJL+Ii- 
1) 1/70. 

AMPL(K)=A0-.25*A1*A1/A2  __  _ 

1400  RETURN 

_ END _ _  ,  .... _ 
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■C10A0  SUBROUTINE  FOR  LOADING  INPUT  DATA  FROM  BUFFER 
SUBROUTINE  LOAD  I  I  SI ZE.EOF ) 

DIMENSION  BUFFER! 12.200.2) .DAT A (12. 207) ,OUT (15 » 200) »B 1 2400.2) 
10(2494) 

COMMON  BUFFER, DATA, OUT. IRUN,STRTSM«FNSHSM»iTRTPR»FNSHPR 
EQUIVALENCE  (BUFFER, B) , (DATA, D) 

_ M*.12*IS  IZE-84 

DO  100  I«l,84 
J«I+M 

100  D( I )«D( J) 

CALL  GET ( JB.EOF , IPFC) 

DO  200  I«X,1PFC 

_ 20CL  D(  1  +  84  )  »B  (  I  ,J9.) 

lS!ZE*7+< IREC/12) 

RETURN 
END 
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3.3  ATTITUDE  DETERMINATION  PROGRAM 


1.3.1  Main  Program.  The  main  program  is  an  executive  program 
which  controls  all  input/output  and  calls  the  main  subprograms  into 
action.  This  program  is  written  in  Fortran  II  for  the  IBM  7094  and 
uses  a  buffered  input/output  on  both  BCD  and  binary  tape  commands. 

The  program  begins  by  initialising  constants  for  the  MAGNET  sub¬ 
routine,  the  detector  angles,  and  associated  trigonometric  functions. 
The  program  continues  and  inputs  a  classified  card,  control  cards  with 
parameters  for  the  specific  rocket,  and  gamma  and  beta  energy  con¬ 
version  tables.  All  inputs  with  the  exception  of  the  classified  in¬ 
formation  are  then  output  for  printing.  The  main  loop  of  the  pro¬ 
gram  starts  at  this  point,  and  the  binary  input  tape  with  smoothed 
magnetometer  data  is  read  into  core  where  the  start  and  stop  times 
are  checked  for  processing. 

Reading  of  the  binary  input  tape  continues  until  the  systems 
time  on  this  tape  equals  the  start  time  indicated  by  the  control  card. 
A  test  is  then  made  to  determine  if  the  systems  time  read-in  falls  in 
the  perturbed  interval  (TLA  to  TLB)  and  is  greater  than  the  atmosphere 
exit  time,  and  switches  are  set  accordingly.  Subroutine  TRAJ  is  now 
called  to  compute  the  coning  axis  vector,  as  a  function  of  the  cur¬ 
rent  systems  time  if  it  is  less  than  or  equal  to  the  atmosphere  exit 
time.  This  procedure  guarantees  that  the  coning  axis  is  constant  above 
the  atmosphere.  TRAJ  is  called  once  more  to  calculate  the  payload 
position  in  Johnston  Island  (JI)  coordinates. 

The  payload  position  is  needed  in  order  to  determine  the  theoret¬ 
ical  magnetic  field  at  that  position  with  the  use  of  the  MAGNET  sub¬ 
routine.  The  MAGNET  subroutine  must  be  given  the  position  in  a 
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latitude-longitude  coordinate  system  and  consequently,  the  first  coor 

dinate  transformation  is  performed.  The  transformation  matrix  is 

defined  in  terms  of  the  launch  latitude  and  launch  longitude  X 

L  L 

east  of  Greenwich  as  follows; 


-sinXj 

cosX. 


-siru|r  cosX 
L  L 

-sinf  sin\ 

L  L 

cost- 


cost^  cosX^ 

cos+  sinX, 
L  L 

siniji. 


tL  °*“tL 

A  vector  x  in  J  I  coordinates  is  expressed  in  geocentric  coor- 
dinates  by 


— / 
x 


G  x 


The  latitude  and  longitude  are  now  computed  using 


%  ‘ Lj,t 


*  lau 


I  '."-j  -  ran 


G-;'"’-? ) 

ffl 


The  elements  of  the  rotation  matrix  which  will  rotate  the  coor¬ 
dinates  of  a  point  in  the  earth's  geocentric  coordinate  system  into 
J  I  coordinates  are  now  computed.  The  resulting  matrix  is  linarized 
since  displacements  are  small  compared  to  one  radian,  and  we  then  ob¬ 
tain  the  following  where  is  the  payload  latitude  and  X^  is  the  pay 
load  longitude, 

1 


B 


(Xp-XL)sintL 

-<Wcos*L 


-<Wsin+P 


p  L 


<Wcos*p 

t  -* 

P  L 


Subroutine  MAGNET  is  now  called,  and  the  components  of  the  theo 
retical  field  x,  in  geocentric  coordinates, are  then  rotated  into  J  I 
coordinates  using  matrix  B,  i.e.. 
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Control  is  now  transferred  to  the  ATUDE  subroutine  where  direction 
cosines  are  determined  and  returned  in  the  form  of  a  matrix  called  A. 
This  matrix  is  now  used  to  compute  the  direction  cosines  of  the  de¬ 
tectors  in  the  payload  which  are  in  turn  used  to  compute  the  detector 
attitudes  using  the  following  equations: 


(Elev.) 


(Azim. ) 


where  Xj ,  y. ,  are  the  x,  y,  and  z  direction  cosines  of  the  ith 
detector,  and  4^  and  <t>i  are  the  elevation  and  azimuth  of  the  ith  de¬ 
tector  . 

The  azimuth  and  elevation  of  the  burst  to  payload  and  rocket 
axis  vector  are  determined  using  similar  equations. 

The  main  loop  of  the  main  program  is  now  completed  by  using  a 
table  look-up  method  of  determining  the  functional  values  of  the 
beta  and  gamma  detectors  from  the  energy  conversion  tables  input  to 
the  program  at  the  beginning. 

All  of  the  computed  results  are  stored  on  a  binary  output  tape,  and 
an  off-line  print  tape  is  prepared  at  this  time. 


GLOSSARY  OF  TERMS 

Fortran  Formula 


Name 

Name 

Description 

FLAT 

+  L 

launch  latitude 

FLONG 

XL 

launch  longitude 

SLAT 

sintL 

CUT 

COS^k 

SLONG 

sin\L 

CLONG 

cosX^ 
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GLOSSARY  OF  TERMS  (contd.) 


Fortran  Formula 


Name 

Name 

DescriDtion 

(XP.YP.ZP) 

-*/ 

X 

geocentric  coordinates  of  a  vector  x 

ELAT 

*p 

latitude  of  payload 

ELONG 

X 

p 

longitude  of  payload 

CONI 

X  -X 

P  L 

factor  of  element  in  lineariz’d  rotation  matrix 

C0N2 

t 

P  L 

factor  of  element  in  linearized  rotation  matrix 

C0N3 

sin* 

P 

factor  of  element  in  linearized  rotation  matrix 

C0N4 

cos* 

„  p 

factor  of  element  in  linearized  rotation  matrix 

FR(I) 

F 

r 

theoretical  field  vector 

DPHI(I) 

*X 

azimuth  of  i  detector 

DTHETA(I) 

«I 

elevation  of  iC^  detector 

A(I ,  J) 

A 

direction  cosine  matrix 

TLA 

start  time  in  systems  time  for  the  use  of  the 
mathematical  model  in  subroutine  ATUDE 

TLB 

stop  time  in  systems  time  for  the  use  of  the 
mathematical  model  in  subroutine  ATUDE 

3-3.2  Subroutine  ATUDE.  Subroutine  ATUDE  determines  the  direction 
cosines  of  the  x,  y,  and  z  magnetometers  (A.j)  in  the  J  X  coordinate 
system. 

This  subroutine  uses  the  coning  axis  vector  CA  and  the  theoretical 
magnetic  field  vector  FR  as  determined  in  the  main  program  and  trans¬ 
mitted  in  common.  The  routine  also  uses  the  subroutine  CONE  to  obtain 
the  coning  parameters  GAMMA  and  BCONE,  the  half-cone  angle  and  total 
azimuthal  coned  angle  of  CA, respectively . 

During  the  magnetic  field  perturbation  period  of  the  Socket 
19  flight,  the  spin  frequency  FREQ  and  spin  phase  angle  PHASE  of  the 
rocket  are  used  in  place  of  the  smoothed  magnetometer  data  in  the  de¬ 
termination  of  the  direction  cosines.  These  two  parameters  are  input 
by  the  main  program  and  transmitted  in  common  to  this  subroutine,  which 
calculates  the  A^  matrix  directly  as  shown  in  the  flow  diagrams  for 
ATUDE . 

The  A^j  matrix  is  transmitted  in  common  to  the  main  program. 
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GLOSSARY  OF  TERMS 


Fortran  Formula 


Name 

Name 

Description 

CA(I),R(I) 

R, 

"3 

coning  axis  vector 

TL 

t 

time  from  launch 

S(I) 

fl 

element  of  e^  orthog.  syst. 

1(1 ) 

e2 

element  of  e^  orthog.  syst. 

BCONE 

0 

the  total  azimuthal  coned  angle  of  R 

GAMMA 

Y 

half  cone  angle 

A(3,I) 

,vv 

,p3,k' 

direction  cosines  of  z  magnetometer 

ATB(I) 

ATC(X) 

“3 

or 

or 

P1 

P2 

^2  for  nonperturbed  period,  P^  otherwise 
U3  for  nonperturbed  period,  ?2  otherwise 

A(1,I) 

P 

A(2,I) 

J1 

FR(I) 

f 

theoretical  field  vector 

FREQ 

rocket  spin  frequency 

PHASE 

6 

rocket  spin  phase  angle 

3.3.3  Subroutine  TRAJ  ■  The  purpose  of  the  TRAJ  subroutine  is  to 
compute  the  position  and  velocity  of  the  vehicle  during  its  flight. 
Since  the  ground  range  of  the  actual  trajectories  is  small  compared 
with  the  radius  of  the  earth,  these  trajectories  are  therefore  approx¬ 
imated  by  parabolas,  i.e.,  gravity  is  assumed  constant  both  in  mag¬ 
nitude  and  direction  over  the  entire  trajectory.  Thus,  the  motion  of 
the  vehicle  can  be  described  in  closed  form.  The  coordinate  system 
employed  is  that  shown  in  Figure  3.1.  The  origin  lies  at  Johnston 
Island,  and  the  x,  y,  z  axes  point  to  the  east,  north  and  vertical, 
respectively.  The  trajectories  are  assumed  to  lie  in  a  vertical  plane, 
containing  the  z  axis,  and  each  is  assumed  to  have  a  constant  azimuth, 
0ff  measured  east  of  north.  The  following  quantities  are  computed,  as 
a  function  of  time,  by  the  subroutine:  x,  y,  z;  x,  y,  z;  x/V,  y/V, 
z/V;  H,  Er>  Here,  V  is  the  velocity  of  the  vehicle,  H  is  the  altitude 

of  the  vehicle  above  the  surface  of  the  (curved)  earth,  and  E  is  the 

r 

earth  range  of  the  vehicle.  These  quantities  are  given  as  follows: 
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where, 


x  -  p  sinflf,  y  -  p  cosa,  z  -  zq  +  At(zQ  -  ^At ) 

X  -  (P0At+(So)sina,  y  -  (poit+(5o)cosa,  z  -  zq  -  gAt 

2  ,  ,  1/2 

V  =  (x  +y  +z  )  ,  E  -  R  tan*1 

r  R+z 

H  .  _ P2+z2+2zR _ 

2  ,  1/2 
Cp  +  (z+r;]  +  r 


2  2 

p  *  (x  4y  ) 


1/2 


Po+At^o42PoAt) 


At 


t  -  t 


R  “  equatorial  radius  of  the  earth  «  6371.2  km 
2 

y  *  9.80  meters/sec 
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The 

begins . 
the  TRAJ 

Fortran 

Name 

ALPHA 

DCVX 

dcvy 

DCVZ 

ER 

G 

H 

R 

RHO 

RHOI 

RHOVI 

RODDT 

T 

TI 

X 

XV 


YV 


time,  tQ,  denotes  the  instant  that  the  trajectory  computation 
Table  3.1  shows  the  value,  of  a,  to,  zq,  zq,  po. Po  employed  in 
subroutine  for  Rockets  8,  9,  15,  19,  and  26. 


GLOSSARY  OF  TERMS 

Formula 

■  Name  Description  and  Physical  Units 

a  angle  (radians)  between  position  vector  and 

Y-axis  (Y-axis  is  assumed  to  point  north) 

cosVx  cosine  (dimensionless)  of  angle  between  velocity 

vectors  and  X-axis 


cosV 

y 

cosV 

z 

earth  range 
g 

H 

R 

P 


cosine  (dimensionless)  of  angle  between  velocity 
vector  and  Y-axis 

cosine  (dimensionless)  of  angle  between  velocity 
vector  and  Z-axis 

range  of  vehicle  (kilometers) 

gravitational  acceleration  (kilometers/second/ 
second) 

altitude  of  vehicle  (kilometers) 

radius  of  the  earth  (kilometers) 

current  value  of  horizontal  position  vector 
(ki lometers ) 


initial  value  of  horizontal  position  vector 
(kilometers) 

initial  value  of  horizontal  velocity  vector 
(kilometers/sec) 

initial  value  of  acceleration  along  horizontal 
(kilometers/ sec^ ) 

current  value  of  time  (seconds) 

initial  value  of  time  (seconds) 

current  value  of  X-coordinate  of  position 
vector  (kilometers) 

current  value  of  X-coordinate  of  velocity 
vector  (kilometers/second) 

current  value  of  Y-coordinate  of  position 
vector  (kilometers) 

current  value  of  Y-coordinate  of  velocity 
vector  (kilometers/second) 


69 


GLOSSARY  OF  TERMS  (contd.) 


Fortran  Formula 

Name  Name 

Z  Z 

ZV  Z 


ZI  Z 

o 

ZVI  Z 

o 


Description  and  Physical  Units 

current  value  of  Z-coordinate  of  position 
vector  (kilometers) 

current  value  of  Z-coordinate  of  velocity 
vectors  (kilometers/second) 

Initial  value  of  Z-coordinate  of  position 
vector  (kilometers) 

initial  value  of  Z-coordinate  of  velocity 
vector  (kilometers/second) 


3.3.4  Subroutine  MAGNET.  The  primary  purpose  of  the  MAGNET  sub¬ 
routine  is  to  generate  theoretical  east,  north, and  vertical  componen 
estimates  of  the  earth's  magnetic  field. 

The  main  field  potential  is  approximated  by  a  truncated  expansion 
of  a  scries  of  associated  Legendre  polynomials 
7  n  n 

v  "  a  Y,  Y  r  pln’n(^)  {F(n,m)  cos(m-l)  0+G(n,m)  sin(m-l)  0} 
n=l  m=l 


where: 


a  =  earth  radius 


r  =  distance  from  earth  center 

0  =  geographic  longitude  (east) 

p  =  cosQ,  0  =  geographic  colatitude 
_  m  *  n  ,  _ 

P  ®  associated  Legendre  polynomials 
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North  Pol* 


\ 

I 


The  gradient  of  V  represents  the  magnetic  field  intensity  vector.  In 
spherical  coordinates  at  the  point  (r,  6,  0)  the  gradient  is  given  as 


W  = 

r  sin8  d0 


dV  -r 
-  1 


.  iav  t  +  sv  ^ 

r  59  J  +  9r 


where  i,  j,and  k  are  the  east,  north, and  vertical  unit  vectors  as 
described  previously. 

7  n  n+l 

BR  =  l7  =  '  1  n  I,  (r)  Pm,n(li)  {F(n,m)  cos(m-l)0 
n=l  m=l 


+  G(n,m) 


B6 


i  av  _  r  Via 
"  r  39  L  L\  r 


n+l 


n=l  m=l 


d(Pm’U(H)  ) 
dp 


|F(n,m)  cos(m-l)0 


+  G(n,m)  sin(m-l)0j 
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7  n 


iltftf  =  Z  Zf  Pn’’n(M')(m-l)  {-F(n,m)  sln(m-l)^ 


n=l  m*l 


+  Gd(n,m)  cos(m-l)$} 


where: 


_in i n .  v  »  >n*l  _ro , n ** 2  _  .  * 

P  00  =  cosQ  P  -  P  C(n,m) 


j  tx^  i  n  .  —I 

-(y)  *  cos0 


,m,  n 


dy 

_ni « n  *  » 
P  (y) 


dp- 


si  n9  P 


i,n-l  dpm'n'2 


dy 


C(n,m) 


/  0,  m  >  n  \ 
i  1 ,  m=n=l  J 


and  C(n,m)  »  (n-2)2  -  (m-1 )2/(2n-3)(2n-5 ) 


C(n,ni)  »  0,  m  >  n 


also 


Pm,n(y )  -  sin0  pra"1,n'1(y) 

dPm,n(y )  A  m-l.n-1 ,  x  „  dPm"1,n‘1(y) 

- — ■ -w-‘-  =  cos0  P  (y)  +  sin©  - - — 

dy  '  dy 


The  constants  F(n,m)  and  G(n,ra)  were  obtained  from  Reference  4 
and  can  be  described  as  arrays  as  follows: 


G(n,m) 


F(n,m) 


0.0 


0.0 

-5798.9 

0.0 

3312.4 

-157.9 

0.0 

1487.0 

-407.5 

21.0 

0.0 

-1182.5 

1000  6 

43.0 

138.5 

0.0 

-79.6 

-200.0 

459.7 

242.1 

-121.8 

0.0 

-575.8 

-873.5 

-340.6 

-11 .8 

-111.6 

0.0 

30411.2 

2403.5 

-3151.8 

-4179.4 

2147.4 

-5125.3 

6213.0 

-4529.8 

-1338." 

-2489.8 

-2179.5 

-649.6 

700.8 

-204.4 

1625.6 

-3440.7 

-1944.7 

-60.8 

277.5 

59.7 

-1952.3 

-485.3 

321.2 

2141.3 

105.1 

22.7 

-32.5 


111  .5 
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The  direction  cosines  of  the  magnetic  vector  then  become 


cos (F,X) 


BR 


2  2  2 
BR  +B0  +B0 


,  cos(F,Y)  =/ 


B9 


2  2  2 
BR+B0+B0 


cos(F,Z) 


B  <)> 


2  2  2 
BR  +B9  +B0 


GLOSSARY 

OF  TERMS - MAGNET  SUBROUTINE 

Fortran 

Formula 

Name 

Name 

DescriDtion 

THETA 

•O 

geographic  colatitude  in  radians 

PHI 

0 

geographic  longitude  in  radians 

R 

Y 

distance  from  earth's  center  in  kilometers 

C 

H,cosO 

S 

sin© 

SP(M) 

sin(m0 ) 

AOR(I) 

(a/r)1 

BR 

BR 

radial  component  of  magnetic  field  (gauss) 

BTHETA 

B© 

northward  component  of  magnetic  field  (gauss) 

BPHI 

B0 

eastward  component  of  magnetic  field  (gauss) 

P(M,N) 

associated  Legendre  polynomials 

DP(M,N) 

dpm,n/dp. 

H 

altitude  in  kilometers 

3.3.5 

Subroutine 

CONE.  This  routine  is  written  in  Fortran  TT  anH 

is  used  for  all  rockets. 

The  routine  determines  for  a  given  time  the  half-cone  angle  GAMMA 
and  the  coning  azimuthal  angle  in  the  coning  axis  coordinate  system 
BCONE . 

For  nonconing  payloads  the  routine  takes  the  form  of  a  dummy  sub¬ 
routine  which  returns  zero  for  BCONE  and  GAMMA. 

For  coning  payloads  an  empirically  determined  coning  buildup 
table  for  GAMMA  and  BCONE  as  a  function  of  time  is  initialized  at  the 
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beginning  of  the  routine.  During  execution  of  the  program,  a  double- 

weighted  parabolic  interpolation  is  used  to  determine  the  values  of  S 

GAMMA  and  BCONE.  If  the  given  time  is  less  than  the  first  entry  in  I 

the  table,  then  this  would  correspond  to  nonconing  time, and  zero  is 

returned  for  both  values  .  | 

During  coning  buildup,  when  the  given  time  falls  in  the  range  of  | 

the  table,  a  single  parabolic  interpolation  is  used  at  the  end  points, 

and  a  double-weighted  parabolic  interpolation  is  used  in  the  subrange  j 

which  consists  of  the  points  between  the  second  and  second-to-last  I 

points  in  the  table.  § 

For  values  of  time  during  uniform  coning,  i.e.,  when  the  given 
time  is  greater  than  the  last  entry  in  the  table,  a  constant  GAMMA  is 
returned  while  the  value  of  BCONE  is  computed  using  a  sawtooth  function. 

f 

This  routine  uses  a  function  subprogram  called  FCN  which  performs 
the  actual  interpolation  using  a  Lagrangian  interpolating  polynomial. 

The  figure  below  is  a  pictorial  description  of  the  general  case  with  the 
corresponding  formulas  where  t  is  the  argument  value  of  time. 


y  (t) 


(Xj-t)  y,(t)  +U-Xj_,)y2(t) 
(Xj  -Xj.,) 
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GLOSSARY  OF  TERMS 


Description 

number  of  points  in  the  table 
systems  time 

used  as  a  switch  and  is  equal  to  one  the  first  time  the 
subroutine  is  called  and  two  every  other  time 

the  independent  variable  in  GAMMA  table  (systems  time) 
the  independent  variable  in  the  BCONE  table  (systems  time) 
the  dependent  variable  in  the  GAMMA  table  (degrees) 
the  dependent  variable  in  the  BCONE  table 
TT 

2n/f;  the  spin  period 

time  at  which  uniform  coning  begins  (systems  time) 
the  half -cone  angle 

the  coning  azimuthal  angle  in  the  coning  axis  coordinate 
system 

3.3.6  Function  ATAN1 ■  This  routine  is  written  in  Fortran  XX  and 
is  a  standard  two-argument  arctangent  routine.  The  subprogram  cal¬ 
culates  the  arctangent  of  y  over  x  and  determines  the  angle  in  radians 
between  -n  and  tt,  i.e., 

-Tt  <  tan  *  ^  £  TT 

x 

The  calling  sequence  is: 

ATAN1  (y  ,x) 

It  is  important  to  note  that  the  arctangent  of  the  first  argument 
divided  by  the  second  argument  is  computed. 

Figure  3.2  illustrates  a  general  case  in  the  second  quadrant  and 
the  angle  9  the  program  computes. 

The  routine  is  used  differently  in  two  places  in  the  main  program. 
The  first  time  it  is  used  to  compute  the  longitude  0  of  the  payload 
using  the  x,  y,  and  z  position  in  geocentric  coordinates  (Figure  3.3). 
This  angle  is  determined  using  the  normal  calling  sequence.  The 


Fortran 

Name 

N 

T 

iswL 

XI  (I) 

X2  (I ) 
Y1(I) 
V2(I) 

PI 

DT 

TU 

GAMMA 

BCONE 
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second  time  the  routine  is  called  it  is  used  to  determine  the  azimuth 
of  the  rocket  axis.  At  this  time  the  calling  sequence  is  reversed  and 
the  routine  computes  ATAN1  (x,y)  (Figure  3.4). 

*i  "  tan'r  7 

By  this  reversal  in  the  calling  sequence  we  are  able  to  compute  both 
angles  with  one  routine. 

3.3.7  Function  FCN.  This  function  subprogram  evaluates  the 
second-degree  Lagrangian  interpolating  polynomial, 

Oc-x2)(x-jc3)  (x-XjXk-Xj)  (x-x^tx-xj) 

f(x'>  (x1-x2)(x1-x3)  yl  +  (x2-Xl)(x2-x3)  y2  +  (x3-x1)(x3-x2)  y3 

which  is  equivalent  to  a  parabola  of  the  form, 
f(x)  =  Ax2  +  Bx  +  C 

through  the  points  {(Xj.y^,  (x2>y2),  (x3,y3)}. 

The  routine  is  written  in  Fortran  II  and  is  used  by  the  CONE  sub¬ 
routine  to  do  the  actual  parabolic  interpolation.  The  calling  sequence 
is 

FCN(xP,x,y,I) 

where  xP  is  the  independent  variable  for  which  the  corresponding 
ordinate  is  desired,  x  and  y  are  arrays  of  the  tabulated  function, 
x  the  independent  variable,  and  y  the  dependent  variable.  The  argu¬ 
ment  I  is  the  subscript  of  the  first  point  to  be  used  from  the  tabu¬ 
lated  function. 

3.3.8  Flow  Diagrams.  The  following  diagrams  outline  the  logical 
and  computational  processes  for  each  of  the  routines  discussed  above 
and  are  followed  by  the  program  coding  as  stated  in  Fortran  and  pre¬ 
ceded  by  a  glossary  of  terms  for  the  convenience  of  the  reader. 
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CALL  TRAJ{T2) 
To  Determine 
Payload  Position 
And  Velocity 
Vectors 


Compute  Detector  Attitudes 


Calculote  Vertical  a  Horizontal 
Gamma  Scanner  Outputs 
In  Counts  /10  m  sec 


Colculoti  Azimuth 
And  El«votion  Of  B*tn 
D»t»c»or,  Horizontal 
a  Vortical  Gamma  Sconnar 
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Compute  Direction  Cotines  Ot 
_Rocket  Axis  Vector  A5 
As=  SGCPS  +  SGSP  T+COSG-R 
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86 


87 


88 


Ent*r  With 


h,  Attitude  in  Kilomtttrt 
i p,  Earth  Latitud#  in  Radiant 
4>,  Earth  Longitud*  in  Radiant 


Stt 

r  =  h+6371.2 


Comput#  Multiple  Angltt 
tin  h<t=tin<&cot(k-l) cj>  -t-cot  <J>  tin(k  — I) 

cot  k$=  cot  4>  cot  (k -l)<X>-tin  <X>tin(k-l) 
K=  2, . 6 


Setup  Initial  Values 
For  BR,  BTHETA,  PHI  Which 
Represent*  The  Three  Components 
After  Accumulations 

BR*  yFO.  U  +  2(*)S  [p<2,  l)TSl  +  P(2,  2)TS2] 
BTHETA*-(i)  [dP(2,  I)  TSI  +  DP(2,  2)  TS2] 

BPHI  *-(-f)3  P(2,  2)  [-P(2,  2)  Sin  ^  +G(2.  2>Co»*] 
a ,  Earth  Rodim 

P«M)*IO  DP(2,I )  *  .  (in.  Q 

P(2,l)«Co«.0  DP(2,2)*  Co*  C 

P(2,2)*  Sin .9 

TSI  •  P(  2, 1 ) 

T  S2  *  F(  2,  2) Cot  <p  +G{2,2)Sin  <f> 
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lit*  Two  Of  Tht  Following 


1.  GAMMA: 

2.  BCONE  : 

3.  GAMMA  = 

4.  BCONE  : 

5.  GAMMA  : 

6.  BCONE  : 


YMI) 

Y2(I) 

FIT, XI, Yl, I) 

F(T,X2,Y2,I) 

(XI  (D-T)'F(T  ,XI,YI,I-2)+(T-XI(I-l))F(T,XI,YI  ,1-1) 
XI(I)-XIU-I) 

(X2(I)-T)F(T,X2,Y2,I-2)  +  (T-X2(I-l»F(T,X2  ,Y2,I-I) 
X2(I)-X2(l- 1  ) 
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AlAhl-f 


Hh  w  ■  jUhi 


GLOSSARY  OF  TERMS 


Fortran  Formula 

Name  Name 

Description 

RAD 

180/tt 

TPI 

2tt 

FLAT 

launch  latitude  (degrees) 

FLONG 

launch  longitude  (degrees) 

F(I,J) 

constants  for  Legendre  polynomials 

G(I,J) 

constants  for  Legendre  polynomials 

THETA(I) 

detector  elevation  angles  (degrees) 

PHI (I) 

detector  azimuth  angles  (degrees) 

STHETA(I) 

sine  of  THETA(I) 

CTHETA(I) 

cosine  of  THETA(I) 

SPHI(I) 

sine  of  PHI (I) 

CPHI(I) 

cosine  of  PHI (I) 

TRIGl(I) 

CTHETA(I)  •  SPHI(I) 

TRIG2(I) 

CTHETA(I)  ■  CPHI(I) 

BIAS 

burst  time  in  systems  time 

XB 

x-coordinate  of  burst  in  J  I 

YB 

y-coordinate  of  burst  in  J  I 

ZB 

z-coordinate  of  burst  in  J  I 

NROCKT 

rocket  number 

TSTART 

processing  start  time 

TSTOP 

processing  stop  time 

GBIAS 

launch  time  in  systems  time 

POTMET 

photometer  multiplying  factor 

TLA 

mathematical  model  start  time  in  L.T 

TLB 

mathematical  model  stop  time  in  L.T. 

FREQ 

spin  frequency 

! 
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GLOSSARY  OF  TERMS  (contd.) 


Fortran 

Formula 

Name 

Name 

Description 

PHASE 

spin  phase  angle 

TA 

atmosphere  exit  time 

ALPHR 

azimuth  of  rocket  trajectory  (degrees) 

GR 

acc  .  of  gravity 

RHOI 

initial  ground  distance 

RHOVI 

initial  ground  rate 

RODDT 

initial  horizontal  rocket  acceleration 

TI 

initial  value  of  time 

ZI 

initial  value  of  z-coordinate  pos .  vector 

ZVI 

initial  value  of  z-coordinate  velocity  vector 

ALPHA 

alpha  in  radians 

SINA 

sine  of  alpha 

COSA 

cosine  of  alpha 

NGAMMA 

number  of  points  in  gamma  table 

VOLTS  (I) 

voltage  table 

GRAMMA(I) 

gamma  energy  table 

NBETA 

number  of  points  in  beta  table 

BNERGY(I) 

beta  energy  table 

SLAT 

sine  of  launch  latitude 

CLAT 

cosine  of  launch  latitude 

SLONG 

sine  of  launch  longitude 

CLONG 

cosine  of  launch  longitude 

IRUN 

run  number  on  binary  X/p  tape 

ICOUNT 

number  of  time  points  read  in 

OUT ( I , J ) 

common  input -output  matrix 

TL 

launch  time 

YM 

y-magnetometer  reading 

XM 

x-magnetometer  reading 

CA(I) 

coning  axis  vector 

ZSUM 

payload  distance  from  center  of  earth 
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GLOSSARY  OF  TERMS  (contd.) 


Fortran  Formula 


Name 

Name 

Description 

XP 

x-coordinate  of  payload  in  geo.  coord. 

syst . 

YP 

y-coordinate  of  payload  in  geo.  coord. 

syst . 

ZP 

z-coordinate  of  payload  in  geo.  coord. 

syst . 

ELAT 

latitude  of  payload  in  geo.  coord,  syst 

. 

ELONG 

longitude  of  payload  in  geo.  coord,  syst. 

CONI 

longitude  of  payload  minus  longitude  of 

launch 

C0N2 

latitude  of  payload  minus  latitude  of  launch 

C0N3 

sine  of  latitude  of  payload 

C0N4 

cosine  of  latitude  of  payload 

H 

altitude  of  payload  from  earth  surface 

BM 

field  magnitude  in  (gauss) 

BR 

radial  component  of  field  (gauss) 

BTHETA 

northward  component  of  field  (gauss) 

BPHI 

eastward  component  of  field  (gauss) 

FR(I ) 

theoretical  field  in  J  I  coordinates 

A(I,J) 

direction  cosines  matrix 

COSDX(I) 

x-direction  of  detectors 

COSDY(I) 

y-direction  of  detectors 

COSDZ(X) 

z-direction  of  detectors 

100 


c 

c 

c 

c 

c 


100 

200 

210 


500 


510 


c 

_ £00 

650 


-SUBROUTINE  ATUQE 

CALCULATES  OIRFCTION  COSINES  OF  X  AND  Y  MAGNETOMETERS 
DIMENSION  A(3,3),CA(3>,CONST(7,7l,F<7,7),FR<3),G(7,7) 

COMMON  A  » ALPHA «C A .CONST  * COSA * F , FR , FREQ » G* GR * ISW1 , ISW2* ISW3.PHASE. 
1RAD.RE.RHOI .RHOVI  .RODDT ,  S  I NA I T  I . TLtXM* YM» 2  1 .2 VI 
DIMENSION  ATB(3)«ATC(3)*P(3).S(3)»TI3) 

-XMODF ( 1 )  *  I-I/A*3 

SWITCH  TWO  . 

GO  TO  (200,100), 1SW2 
TURN  SWITCH  THREE  ON 
I SW3  «  1 

COMPUTE  COMPONENTS  OF  E 1  BAR  VECTOR  - - 

DO. 210  1*1,3 

R(  I)  =  CA(  I)  '  - - 

SWITCH  THPE.E 

GO  TO  (500,900 ) , I SW3 

COMPUTE  DOT  PRODUCT  OF  MAGNETIC  FIELD  AND  VELOCITY  VECTORS 
DOTFV  =  0.  “■ 

-00.  5.10_  1*1,3 

DOTFV  =  DOTFV  +  FR(I)*CA(!) 

COMPUTE  SINE  OF  ANGLE  BETWEEN  THESE  VECTORS  AND  CHECK  MAGNITUDE 

SIN2  *  SORTF ( 1 , -DOTFV** 2 )  - 

I F(SIN2-l.E-5 >600,600, 700 

ERROR  . . 

PRINT  650, Tl 


JL 

C 

~C~ 

C 

c 


800 

900 


^FORMAT(///10H  WHEN  TL  *E1«.7,34H,  CONING  AXIS  COINCIDES  WITH  FIELD 

GO  TO  1700  - 

700  DO  710  1*1,3 

710  S(I)  =  <FP( I )-DOTFV*R(  I  > ) /SIN2  ~ 

COMPUTE  COMPONENTS  OF  E2BAP  VECTOR 

DO  720  1  =  1,3  . . 

J  *  XMODFU  +  1) 

K  *  XMODF (1+2)  - 

720  T ( I )  *  R(J)*S(K!-R(K)*S(J) 

SWITCH  TWO  - - 

GO  TO  (800,900) , ISW2 

TURN  SWITCH  THREE  OFF  "  - 

I SW3  =  2 

CALL  CONFITL, BCONF, GAMMA, ISW1  )  . . 

TURN  SWITCH  ONF  OFF 

I SW1  =  2  - - 

BEGIN  ACTUAL  COMPUTATIONS 

GAMMA  =  GAMMA/57.29578  - - - 

CO SP  *  COSF(BCONE) 

SINP  *  S  INF  (BCONE  )  - - 

COSG  =  COSF(GAMMA) 

SING  «  SINF(GAMMA) 

SGCP  *  S I NG*COSP 

SGSP  =  SING*SINP  - - - 

COMPUTE  DIRECTION'  COSINES  OF  ROCKET  AXIS  VECTOR 
DO  O10  1*1,3 

A (3, I)  =  SGCP*S( I ) +SGSP*T ( I )+COSG*R ( I ) 

SWITCH  TWO 

GO  TO  (1500,1000) ,ISW2 

COMPUTE  DOT  PRODUCT  OF  ROCKET  AXIS  AND  MAGNETIC  FIELD  VECTORS 
’*#2+A<  3*2>  **2+A(  3»3>**2-,9  11025, 1025, 1050 

PRINT  1026, TL 


910 

C 

c"' 

1000 

C1025 

-Ifl.25 
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_J026  [!°R”£T^0WHEN  TU?E14.7«34H,  ROCKET  AXIS  IS  NOT  A  UNIT  VECTOR) 

1050  AFOOT  *  0. 

00  1050  1=1,3 

1060  AFOOT  =  AFPOT+A  (  3, 1 ) *FR( I ) 

SIN1  =  SORTF ( 1 , 0- AFOOT* AFOOT ) 

IF(SIN1-1.F-5)1  TOO, 1200, 1?00 

C  FRROR  '  - 

1100  PRINT  1150, TL 

1»°  JOJNAT,„/10H  WHEN  TL  =E14.7,35H,  PAYLOAD  AXIS  COINCIDES  WITH  FI  EL 

GO  TO  1700  -  . 

C  COMPUTE  COMPONENTS  OF  U2BAR  VECTOR 

1200  00  1210  1  =  1  ,3  - 

1210  ATBII)  =  ( FR ( I ) -AFOOT* A ( 3 , 1  1  I /S I N1 

C  COMPUTE  COMPONENTS  OF  U3BAR  VECTOR  - 

00  1 770  1  =  1  ,3 

J  «  XMOOFII+1)  -  '  — ■ 

_ K  =  XMODF (1+7) 

1220  ATC(I)  =  A(3,J)*ATR(K)-A<3*K)*ATB(J)  ' 

C  COMPUTE  MAGNETIC  FIELD  VECTOR  IN  X-Y  PLANE 

SMTR  =  SORTF  (XM#*2+YM*«2  )  . . 

IFISMTR-l.E-5) 1300, 1400,1400 
C  ERROR 

_ 1.300  PR  I  NT  1350.TL 

1350  FORMAT ( 9H  WHEN  TL=E 1 4. 7 , 23H ,  SMTR  IS  LESS  THAN  1-51  — 

GO  TO  1700 

1400  COSCHI  =  XM/SMTR  . . 

S INCHI  =  -YM/SMTR 

GO  TO  1600  -  '  - 

_L5_Q0_CGCP,=  COSG*COSP 

CGSP  =  COSG*S I NP  . 

00  1510  1=1  ,3 

ATB(I)  =-S ( I  I  *S I NP+  T I I ) *COSP  - - 

1510  ATC(I)  =  SING*R( I )-(CGCP*S( I )+CGSP*T(I) I 

ARG  =  FREO*TL-PHASE  '  - 

_  COSCHI  =  COSF(ARG) 

SINCHI  =  SINFIARG)  '  '  - 

C  COMPUTE  DIRECTION  COSINES  OF  X,Y  MAGNETOMETERS 

1600  00  1610  1  =  1  ,3  - — . - 

All, I)  =  C0SCH1#ATRI I)+SINCHI*ATC( I) 

1610  A 1 2 , 1 )  =  COSCHI*ATC(I)-SINCHI»ATB(I)  '  . ~ 

GO  TO  IS 00 

C  SET  DIRECTION  COSINES  TO  ZERO 
1700  00  171 0  1  =  1  ,3 
DO  1710  J=1  ,3 
1710  A  I  I , J )  =  0. 

1800  CONTINUE 
RETURN 

END  ■  - 
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.SUBROUTINE  TRAJU».X±.Y,ZtDCVX*DCVY*DCVZ.ER,H)  . . 

DIMENSION  A (3,3) ,CA< 3 ) .CONST ( 7,7 ) ,F(7,7),FR(3),G(7.7> 

COMMON  A, ALPHA, CA, CONST,  COSA,F,FR,FREQ»G,GR,ISWl • I SW2 . 1 SW3 .PHASE, 
1RAO.RE.RHOI .RHOVI.RODDT.SINA.TI ,TL»XM»YM.ZI ,ZVI 
VMTI  =  T-TI 

DETERMINE  MAGNITUDE  OF  HORIZONTAL  POSITION  VECTOR 

RHO. .*  RHOI  +  TMT I*(RHOVI  +  » 5*RODDT*TMT 1 1  _  _ _ _ 

DETERMINE  X.*,  AND  7  COORDINATES  OF  POSITION  VECTOR 
X  »  RHO*S!NA 

Y  «  RHO*COCA 

Z  ■  Z I+TMT I* ( ZVI-,  5*GR*TMT  I  )  .  ...... 

DETERMINE  XV, YV,  AND  ZV  COORDINATES  OF  VELOCITY  VECTOR 

XV.  ■  (RODDT*TMTI  +  RHOV I ) *S I NA  _ 

YV  »  COSA*( RODDT *TMT I  +  RHOVI ) 

ZV  «  Z VI -GR*TMT I 

DETERMINE  MAGNITUDE  OF  VELOCITY  VECTOR 

V  »  SQRTF(XV**2+YV**2+ZV**2) 

DETERMINE  DIRECTION  COSINES  OF  VELOCITY  VECTOR 

DCVX  «  XV/V  _  _ 

DCVY  -  YV/V 
DCVZ  =  ZV/V 
DETERMINE  ALTITUDF 

H  *  (RHO**2+Z**?+2.*Z*RE)  '  ( SQRTF ( RHO**2+ ( Z+RE > **2 ) +RE ) 

DETERMINE  FARTH  RANGE 

ER  *_  RE* 1  TANF(RHO/(7  +  RE')  _ 

RETURN 

END 


nnnnn 


SUBROUTINE  MAGNET ( H, ELAT , £LONG,B,BR* BTHETA  *BPH I ) 

DIMENSION  A (3. 3) , CA ( 3 ) , CONS T ( 7 ,7  I , F ( 7, 7  I ,FR 1 3  1 ,61 7 ,7 ) 

COMMON  A* ALPHA, CA  * CONST , COS A. F ,FR, FREQ* G.GR* ISW1 , ISW2 , I $W3 , PHASE . 
1RAD,RE,RH0I »PHOVI • PODDT ,S I NA  ,  T I ,TL,XM,YM*ZI*ZVI  >5wa. PHASE. 

DIMENSION  AOP ( 8 ) , CPI 7 ) ,DP 1 7,7 1 ,P< 7,7 » • SP ( 7 ) 

C  H  ■=  ALTITUDE  IN  KILOMETERS 

_C -  ELAT  «  EARTH  LATITUDE  IN  RADIANS 

ELONG  *  EARTH  LONGITUDE  IN  RADIANS  .  . 

B  =  MAGNETIC  FIELD  MAGNITUDE  IN  GAUSS 
BR  =  RADIAL  COMPONENT  IN  GAUSS 
BTHETA  *  NORTHWARD  COMPONENT  IN  GAUSS 
BPHI  =  EASTWARD  COMPONENT  IN  GAUSS 
_  THETA  =  1. 57079634-FLAT 

PHI  »  ELONG  - 

R  *  H+6371.2 
C  *  COSE (THETA  I 
S  =  SINF(THFTA) 

C  MULTIPLE  ANGLE  FORMULAE  . 

_ SPM  1=0,0  .. 

CPin-i.o  . — ■ — - 


SP< 2 )=SINF(PHI) 

CP ( 2 ) =  COSFIPHI  | 

DO  50  M=3,7 

SP(M)  =  SP ( 2 >*CP (M-l )+CP( 2 )*SP(M-1  ) 

- .50  CP  ( M )  =  CP(2)*CP(M-n-SP(2)*SP(M-l) 

C  POWERS  OF  A/P 


AOR(I)  *  6371. 2/R 
DO  60  I  *  2.8 

60  AOR(I)  *  AOR (1 )*AOR< 1-1  I 
BR  =  0. 

-  BTHETA  =  0. 

BPHI  =  0. 

DO  70  L-»1.7 
DO  70  l 2=1,7 
PIL1.L2I  =  0. 

70  DPI  LI ,L2 1  =  0. 

_ pu.i.i  =  1.0 

PI2.1  )  =  C 
P  1 2 . 2  »  =  S 
DPI  2 , 1 >=-S 
DP ( 2 , 2  )  *C 
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90 


TS1=  F ( 2 , 1 1 *CP ( 1 ) +  G I  2, 1 ) #SP  1 1 ) 

TS2=  F (2,2 )*CP( 2 >  +  G(2,2)#SP(2t 

BR=  AOR (2 ) »F ( 1 ,1  )  +  2.0*A0R(3)*(P(2,l )*TS1+P(2,2)*TS2 
8THETA=-A0R ( 3 1  * ( DP ( 2,1 )*TS1+  DPI2,2)*TS2I 

BPHI=-AOR(3>*P(2,2)*(-F(2,2)#SP(2)+G(2,2)*CP(2l ) 

on  i,n  n  =  3,7  . 


DO  120 
FN  =  N 
SUMR  =  0. 
SUMT  * 
SUMP  s 

do  no 


0. 

o, 

M 


=  1  ,N 


GENERATION  OF  ASSOCIATED  LEGENDRE  POLYNOMIALS 
IE  (N-M)  90,80,90 
K  =  M-l 

PIN, Ml  »  S*P (N-l ,K  J 

DP ( N , M  )  =  S*DP (N-1,K)+C*P( N— 1 ,K I 

GO  TO  100 

P(N,M)  =  C»P (N-l ,M) -CONST (N,M) «PIN-Z,M I 
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_ _  OP(N,M)  •  C*0P(N-ltM)-S*P(N-l,M)-CONST<N,MI*DP(N-2.M) 

100  CONTINUF 
FM  •  M-l 

C  COMPUTE  SUMMATIONS 

TS  «  F ( N |M) *CP ( M ) +G I N  *M )*SP ( M I 
SUMR  «  SUMP+P (N«M)*TS 

_ SUMT  «  SUMT+DP ( N  *M ) *TS 

liO  SUMP  «  SUMP+FM*P(N.M)*(-F(N.MI*SP(M)+G(N*M)*CP(MM 

BR  »  BR+AOR (N+l  I  *FN*SUMP 
BTHETA  *  BTHETA-AOR I N+l )*SUMT 
BPH I  «  RPHI-AORfN+l ) *SUMP 
120  CONTINUE 

_ BPH  I  ■  BPHI/S  ..  _ 

B  ■  SQPTF ( RP*BR+STHETA*BTHET A+BPHl *BPH I ) 

C  SWITCH  FROM  JENSEN-CAIN  SYSTEM  TO  EOS  SYSTEM 

BR  *  -RP 
BPHI  *  -BPHI 
RFTURN 
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XCQNE  ..SUBROUTINE  CONE  FOR  ROCKET  NO.  19 
SUBROUTINE  CONE  (TL.BCONE. GAMMA. ISW1 I 
DIMENSION  X1(7),X2(7)»Y1(7),Y2(7I 
N*7 

T  *  TL+7686.4 
GO  TO  (10,30) ,ISW1 
_  _  10  XI ( 1 1*7735. 

Xl(2)*7739.72 
X 1 ( 3 1*7745.57 
Xl(4)«7752.17 
X 1 ( 5 1 *7761 .91 
XI (61*7774.59 

___ _ Xl<7)*7787.29 

X2< 11*7737.44 
X2 ( 2 1 *7742 , 

X2(3)*7749.14 
X2(4)*7755.21 
X2(5)*7768.61 
____X2  (61*7780 ._56 
X2 ( 71*7794,01 
Y1  ( 1 1 *0, 

Y1  <  2 1 *3 .4 
Y1  ( 31*4.916 
Y1  (41*2.90 

_ _ X115 1  *.?,!  1.7..... 

Y 1  (61=7.634 
Yl(7)*7,936 
P I  =  3. 1  41 5926R 
DT  *  2.*PI/.2439 
TU  *  X  2 ( N 1 

_ DO  20  1*1, N 

XI-I-1 

20  Y2(I)*XI*PI 
30  IFIT-Xl (1  1 140,40,50 
40  GAMMA=0. 

GO  TO  65 

_ SO..  J E (  T.-X 2 ( 1 1)60,60,70 

60  GAMMA*FCN  (T.Xl.Yl.l) 

65  BCONF*0. 

GO  TO  310 

70  I F ( T-Xl (21 180,90,100 
80  GAMMA»FCN  (T,X1 ,Y1  ,1  1 

_  GO  TO  95 

90  GAMMA=Y1<2) 

95  BC0NE*FCN  (T,X2,Y2,1) 

GO  TO  310 

100  IFIT-X 2(2)1110, 120, 130 
110  BCONE=FCN  (T,X2,Y2,1) 

.  ....  GO.  TO  125 
120  BCONF=Y2 ( 2 1 

'l25l-XM2M  <X1 ,31~T,*FCN  ,T’X1’Y1,1,  +  ,T“X1 ,2> ,#FCN  (T,X1,Y1»2)>/(X1(3) 
GO  TO  310 
130  M=N-2 

00  140  1=3, M 
IFIT-Xl  ( 1 1 1260,270,135 
135  !F(T-X2( I  1 1280,290,140 
140  CONTINUE 
1=N-1 
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1FIT-X1.<II)260, 160*170 
160  GAMMA=Y1 ( I ) 

GO  TO  285 

170  I F  t  T-X2  1 1 1 )  180* 1 90* 200 
180  GAMMA-FCN  (T.Xl.Yl.M) 

GO  TO  285 
—  190  BC0NE=Y2(I) 

GO  TO  215 

200  1F(T-X1 (Nl  1210,220.230 
210  BCONE'FCN  (T,X2.Y2,M> 

215  GAMMA=FCN  IT, XI, Yl, Ml 
GO  TO  310 

_ 220,  GAMMA* Y.1  (N) 

GO  TO  245 

230  IF(T-X?CN> 1240,250,300 

240  GAMMA=Y1 ( N I 

245  BCONE*FCN  (T,X2.Y2,M) 

GO  TO  310 

_JL50_.BCQNE*Y2.(„N)  .... _ 

GAMMA  *  YKNI 
GO  TO  310 

260  GAMMA* ( CXI ( 1 )-T )*FCN  { T »X1 , Y1 , I -2 >♦ ( T-Xl II -1 ) ) *FCN 

i<xi<n-xni-m 

GO  TO  285 
__270  5AMMA  =  YUn 
GO  TO  285 

280  GAMMA* ( < X 1 (  1  +  1 ) -T ) *FCN  ( T ,X1 , Y1 *  1-1)  +  ( T-Xl ( I ) ) *FCN 
11 ( 1  +  1 J  —X 1  (  1 1 ) 

285  BCONE* t ( X2< I >-T I *FCN  ( T .X2.Y2 , I -2 1  + ( T-X2 ( 1-1 )) *FCN 
l(X2(n-X2(!-ln 

. . ..GO  TO  310  .  ..  . 

290  BC0NE*Y2  1 1  1 

GAMMA* ( <X1(  1+1 )-T)*FCN  ( T ,X1 , Y1 ,  1-1 )  + 1 T-Xl ( It )*FCN 
11 ( 1+1 1 —X 1 ( I) 1 
GO  TO  310 
300  GAMMA*Y1(N) 

_ 302  I FJLT-J  TU+DT  1 1306*306,304 

304  TU  =  TU  +  DT 

306  BCONE  *  • 2439» ( T -TU) 

310  RETURN 
END 


(T*X1,Y1, 1-1! )/ 


CT,X1,Y1,XU/IX_ 

(T$XZfY2tJ-lJU- 


(T»X1,Y1»JJ  USJL 
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FUNCTION  ATAN1(Y,X) 

C  THIS  S  U  B  R  0  Uf IN  ECA  L  CU  LA?  ET  THE  A  ft  CTA  NG  ENT-  OF  (Y/X»  IN  THE  INTERNAL 

c -PI  TO  +PI.  IF  V  ANO  X  ARE  ZERO,  THEN  ATAN1F  ( Y  »X )  IS  SET  TO  ZERO. 

IF (Y) 10,70,10'  .  . - 

10  IF (X ) 20 , 30 , AO 

20  AT AN  1*3.14159 26 8-ATANFI  ABSFTyTxTT 

_ Sfl.  TQ  5_Q _ 

30  ATAN1*1. 57079634 


40  ATAN1*ATANF  (ABSF  ( Y/X  )  ) 
60  ATAN 1*-ATAN 1 


_ Qo_mioo _ 

70  IF  (X  180 ,90 , 90 

_ &0_  AI  AH1= 3  ..1 41 5  9  268 . 

GO  TO  100 

_ 90  AT  AN  1*0. _ _ 

100  RETURN 

_ END _ 
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FUNCTION  FCN(XP,XOYT+) 

DIMENSION  XTTTTY ( H 
XHOOF(I)  =  I-I/4*3 

— “rarw; - - — — — 

OO  20  J-1,3 

H,I+J.1  - 

YP*Y (M  ) 

DO  10  1*1,2  - - 

N  -  I+XMODF ( J+L )  -1 

— rO'YP«YPV(XT>-X  XTW  T -XTNTT - 

20  FCN-FCN+YP 
RTTCfftN 
ENO 
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CHAIN 


C 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


ATTITUDE  DETERMINATION  PROGRAM 

DIMENSION  A(3,3),CA(3),CONST(7,7),F(7,7)  ,FR  ( 3 )  (  7 , 7  )  ,QUT  146,200  ) 

COMMON  A, ALPHA, C A,  CONST, COS A, F,  FR, FREO,G,GR, ISW1 ,  ISW2, I SW 3, PHASE, 

1RAD,RE,RHOI,RHOVI,RODDT,SINA,TI,TL,XM, ym.zj , zv I . 

DIMENSION  BNERGYI68 ) , CDSDXI 3 ) ,  COSDY ( 3 ) , COSDZ ( 3 ) , C  PH  I ( 3 ) , CTHE  T A ( 3) 
1DPHI (3) .DTHETAi 3) , GR AMMA ( 6 6 ) , PH  1 1  3 ) , SCOSQI 3 ) ,SPHI ( 3 )  ,  STHET A ( 3 ) , 
2THET  A I  3 ) , T  R I G 1 ( 31 ,1RIG2(3),VULTS<68) 

TSTART  =  TIME  TO  START  PRUCESSING  IN  SECONDS  FROM  LAUNCH 
TSTOP  =  TIME  TO  STOP  PROCESSING  IN  SECONDS  FROM  LAUNCH 
BIAS  =  BURST  TIME  IN  SYSTEM  TIME 
GBIAS  =  LAUNCH  TIME  IN  SYSTEM  TIME 
FLAT  =  LAUNCH  LATITUDE  IN  DEGRtES 

FLONG  =  LAUNCH  LONGITUDE  IN  DEGREESIEAST  UE  GREENWICH) 

BAAL  =  BURST  ALTITUDE  IN  KILOMETERS 

N GAMMA  =  LENGTH  OF  GAMMA/VOLTS  TABLE 

NBETA  =  LENGTH  OF  BETA  ENERGY  CONVERSION  TABLE... 

PUTMET  =  PHOTOMETER  MULTIPLYING  FACTOR 

BCONE  =  TOTAL  CONING  ANGLE  RADIANS 

FREQ  =  SPIN  FREQUENCY  IN  RADIANS  PER  SECOND 


GAMMA  =  HALF  CONE  ANGLE  IN  DEGREES 
PHASE  =  DELTAIPHASE  ANGLE ) I N  RADIANS 
THEi'A(l)  =  ELEVATIuN  ANGLE  OF  VGS20  IN  DEGREES 

THE T A (  2 )  =  ELEVATIUN  ANGLE  OF  HGS20  IN  DEGREES 

THET  A  (  3 )  =  ELEVATIUN  ANGLE  OF  BETA  DET  ECTOR/ PHOTOMETER  -  IN  DEGREES 

PH  1(1)  =  AZIMUTH  ANGLE! FROM  Y)  OF  VGS20  IN  DEGREES 
PHI  (2)  =  AZIMUTH  ANGLE  OF  HGS20  IN  DEGREES 

PHI  (3)  «  AZIMUTH  AinGLE  OF  BETA  DETECT  OR  /  PHOTDMET  E  R  IN  DEGREES 
SWITCH  DISCRETION 
SWITCH  ONE 


ON  IF  SUBRUOTINE  CONE  HAS  NOT  BEEN  CALLED 
OFF  IF  SUBROUTINE  CONE  HAS  SEEN  CALLED 
SWITCH  TWO 

ON  IF  TL  IS  IN  (TLA, TLB) 

OFF  IF  .11..  IS  NOT  IN  (TLA, TLB) 

SWITCH  THREE 

ON  IF  RI..AND  FRI  VECTORS  HAVE  CHANGEO 
OFF  IF  R I  AND  FRI  VECTORS  ARE  CUNSTANT 
9  FORMAT (BE9. 2) 

10  FORMAT  (6612.8) 

11  FORMAT (I2/I8E9.2)) 

12  FORMAT ( I  2/ (  16F4. 2  )  ) 

13  FORMAT  (I2/I6E12.8) ) 

SET  UP  CONSTANTS 
RAO  =  57.2957800 

TP  I  =  6.28318536 
FLAT  =  16.734250 
FLONG  =  169.528189 
NPAGE  =  0 

RE  =  EARTH  RADIUS  IN  KILOMETERS 
RE  =  6371.2 
F(  I,  1  I  =  0. 

F  (  2 ,  1 )  =  30411.2 
F ( 3 , 1 )  =  2403.5 
F  (4, 1)  =  -3151.8 
F  (  5 , 1 )  =  -4179.4 
F  ( 6 , 1 )  =  1625.6 
F(7,  I )  =  -1952.3 
F ( 2, 2 )  =  2147.4 
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F  ( 3 , 2  ). 

= 

-5125.3 

F  (4 , 2 ) 

= 

6213.0 

F  ( 5 , 2 ) 

= 

-4529.8 

-  ■ 

F  (6 , 2 ) 

= 

-3440.7 

F .( 7 , 2 ) 

= 

-485.3.  .  . 

. .  . . .  — - 

- . .  -- 

F  {  3 , 3  ) 

= 

-1338.1 

F  (4 ,3  ) 

S 

-2489.8 

. -  ■ 

- . - . 

FIS, 3) 

= 

-2179.5 

F<6,3> 

= 

-1944.7 

“* 

. 

F  ( 7 , 3  ) 

= 

321.2 

FI4.4  ) 

5* 

-649.6  - 

-  . — . 

FIS, 4) 

= 

7  00.8 

F  (6 ,4  ) 

= 

-60.8 

. -  • 

-  —  • 

F I  7 ,4  ) 

= 

2141.3 

FIS, 5) 

s 

-204.4 

- 

•  - 

F  1 6 , 5  ) 

= 

277.5 

F  17 ,5  ) 

= 

105.1  . . 

-  . . 

- - 

F 16,6  ) 

a 

69.7 

F  1 7 ,6  ) 

8 

22.7 

•  -  -  ~  • 

. 

F  1 7 ,7  ) 

= 

111.5 

G I  1,1) 

= 

0. 

... 

G I  2,  1 ) 

= 

0. 

G 1 3, 1 ) 

s 

o. . 

— 

GI4, 1) 

S 

0. 

G 1 5  *1  > 

s 

0. 

•  •  -  - 

. 

G  1 6 , 1 ) 

s 

0. 

G 1 7 , 1 ) 

s 

0. 

. .  *-  * 

G  I  2 , 2 ) 

B 

-5798.9 

G 13.,  2 ) 

B. 

-  33 12  .4 - - -  -  - 

- -  - . 

GI4.2) 

= 

1487.0 

G15.27 

B 

— 1182.5- 

....  ..  - - - 

. . . 

G  (6 , 2 ) 

3 

-79.6 

G 1 7 , 2  ) 

3 

-575.8 

... 

"  1  --  ‘  ' 

G  1 3 , 3  ) 

B 

-157.9 

G(4  *3  ) 

.  B 

-407.5 

•  —  -  - - - . -  •  • 

...  —  ... 

G(5,3) 

= 

1000.6 

GI6.3  1 

= 

-200.0  . 

- - - - 

—  — 

G  1 7 , 3  ) 

= 

-873.5 

G 1 4 ,4  ) 

= 

21.0 

G 1 5 ,4 } 

= 

43.0 

G 1 6 ,4 ) 

= 

459.7 

— - . —  -  - 

....  ..  -  - 

G 1 7 ,4 ) 

3 

-340.6 

GI5.5) 

= 

138.5 

-  -  —  •  • 

. 

G 1 6 , 5  1 

a 

242.1 

G 1 7 , 5  ) 

B 

-11.8 

.. 

.  . 

G  16,6  ) 

s 

-121.8 

‘,17,6  ) 

= 

-111.6 

...  . . - . -  - 

-  -  —  ■  - 

G  1 7  ,7) 

~ 

-32.5 

00  20 

1  = 

1,7 

-  - 

-  ..  - 

DO  20 

J= 

1,1 

FI  I , J) 

s 

F 1 1 , J ) *  1 .OE-5 

... 

- 

20  G 1 1 , J ) 

= 

G  ( I  ,  J  )  *1. 0E-5 

C  COMPUTE  CONST  .  .  .  . 

00  30  N  =  1,7 

FN  *  N  -  - 

00  30  M  =  1,N 

FM  =  M  ..... 

30  CONST  (N,M  )  =  (  I  FN-2.0  I  **2- (  FM-1 .0  ) #*2 )  / (  (  FN+FN-3.0 )  *<  FN  +  FN-5 .0  ) ) 
C  INITIALIZE  DETECTOR  ANGLES  -  ■  - 
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THETA!  1 )  =  0.0 
THE T A ( 2  )  =  -15.0 
THETA! 3 )  =  0.0 
PHI  (  1 )  =  0.0 
P_H I  (  2 )  =  45.0... 

PHI (3)  =  90.0 
DO  15  I  =  1,3 

STHETA(I)  =  SINF  I  THETA!  I  )  /RAO) 
CTHETA(I)  =  C0SF(THETA(1  l/RAO) 
SPHI  II)  =  SINF (PH  I  I  I (/HAD) 

C.PHI  (I  )  =  COSFIPHI  (I  )  /RAO) 
TRIGKI  )  =  CTHETA  ( I  )*SPHI  (I  ) 

15  TRIG21I)  =  CTHETA! I)*CPHI < 1) 

DO  17  1=41,68 


ENERGY ( I ) =0. 

VOL  T  S ( I )=0. 

17  GRAMMA! 11=0. 

TURN  SWITCHES  ONE  AND  THREE  ON 
ISW1  =  1 

ISW3  =  1  '  '  “  . - 

READ  CONTROL  CARD 

READ  INPUT  TAPE  5  , 10 , 81 A  S ,  XB ,  YB ,  ZB 

INPUT  ROCKET  NUMBER.  AND  ASSOCIATED  CONSTANTS. _  _ 

RE  AO  INPUT  TAPE  5 , 1 3,  NROCKT ,  TSTART ,  TSTOP ,  GB I  a’s  TpOTMET  ,  T  L  A  j  TL  B~,  F  RE  0 
1 , PH ASE , TA 

INPUT  CONDITIONS  FOR  ROCKET  TRAJECTORY 

READ  INPUT  TAPE  5,  10  ,  ALPHR  ,  GR  ,R  HO  I  ,  RHOVI  .RODDT  ,J  I  ,  Z 1  ,  Z.V  I  . 

CONVERT  ALPHA  TO  RADIANS  AMD  COMPUTE  SINA  AND  COSA 
ALPHA  =  ALPHR/RAD  .  .. .  _ 

SINA  =  SINF!  ALPHA  ) 

COSA  .=  COSF  (ALPHA) 

READ  IN  TOTAL  GAMMA  DETECTOR  ENERGY  CONVERSION  TABLE . 

READ  INPUT  TAPE  .5,12,  NGAMMA,  (VOLTS!  I),  1  =  1,  NGAMMAJ 
READ  INPUT  TAPE  5  ,9 ,  ( GR  AMMA  ( I  ) ,  I  *  1 ,  NGAMMA  ) 

READ  IN  BETA  CONVERSION  TABLE 

READ  INPUT  TAPE  5  ,  1 1 ,  NBE  T  A  ,  (  BNERG  Y  (  I  )  ,  I  =1  ,  NBE  TA  ) 

OUTPUT  INPUT  LIST 

WRITE  OUTPUT  TAPE  6 , 35 , MROCK T, TSTART, TLA , POTMET, TS TOP , T L B  , FREO , GB I 
IAS ,T  A, PHASE 

35  FORMAT  (  1H  144X28HINPUT  DATA  FOR  ROCKET  NUMBER  I  3//  1 H08X  1 BHCONTROL  PA 
1REMETERS/1H013X8HTSTART  =F5.1,4H  SEC13X5HTLA  =  1PE8.1,4H  SEC  13X8HPO 
2TMET  =E9. 2,  26 H  W A TT /CM** 2/ S TERAD I  AN/VOLT/ 1  HOI  3X7 HTSTOP " =0PF6. 1 , 4H 
3SEC13X5HTLB  =1PE8.1,4H  SEC13X6HFREQ  =E14.8,8H  RAD/S  EC /JHO 1  3X7HGB I A 
4S  =E11.4,4H  SEC8X4HTA  =0PF10.4,4H  S EC 12X7HPHA SE  =1PE14.8,4H  RAD//) 

WRITE  OUTPUT  TAPE  6,  36,  AL  PHR ,  T1 ,  RHO I ,  RODDT  ,  f.LATjRE  ,  RHOY  I  ,GR  i  F  LONG , 

1 Z  I  t  Z  V  I 

36  FORMAT  (  1HG8X  32HC0NDI  TI  ONS  FOR  ROCKET  T  RA  JEC  I'ORY/ 1  HO  13X7H  AL  PHA  =F  10 

1.5 , 4H  DEG 10X4HTI  =F6.1,4H  S EC  1 1 X7HRH0 I  =F7.3,3H  KM14X7HR0DDT  =1PE 
2U.4/1H013X7HFLAT  =0PF10.5,4H  DEG10X4HRE  =F7.1_,3H  KM  1  J  X  7HRH0V  I  = 

3F7 .4 ,7H  KM/SEC10X7HG  = 1  PE  1 1 . 4/ 1  HOI 3X7HFLONG  =0PF10.5,4H  DEG10X 
44HZI  =F8.3,3H  KM10X7HZVI  =F7„4,7H  KM/SEC//7 ) 

WRITE  OUTPUT  TAPE  6,37  .  . 

37  FORMAT (49X 21HBETA  CONVERSION  TABLE/48X23H! FLUX  IN  M EV /CM** 2/SE C ) // 
14X6 ! 5HV0LTS4X4HFLUX6X ) // ) 

DIMENSION  P  (  6  ) 

DO  38  1=1,6 

38  P ( I ) =0 . 

DO  40  1  =  1, NBE  TA,  6  .  ~  ■  -  -  ■ 

P( 1)=P(6)+. 1 
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Q 

C 


.53 

C 

..55 


39  P(J)«P( J-l)+.l  -  -  . 

p?J>ULIor!l?,^,P(U,BNERGY'  1,1  P(  21  -BNtRGY  I  1  + 1 )  ,P(3>,  BNERGY 
.4 INFORM AT  ^MOPF  9^  B*'*ER6  I+/* I  » P(  6  )  ,BNERGY  (  I  +B  ) 

WRITE  OUTPUT  TAPE  6,42"  "  '  *  '  '  . -  —  -  - 

■■■42ic?mIif5^?i«w??uS2«55fION- ^TABLE/48x23H<f‘-UX.INj)6yyW**2/SE 
iSn/Jc  ,?UI?UT  TAPE  b,41‘  1  VOL  TS  { I  ),  GRAMMA  t  I)  .VOLTS  (1+1)  , GRAMMA  (  1  +  1). 

. i4?tvlLT;!i;i^GRAMA:fi;5?^=i!NliKSA^r,,+3,’!_0LTS<I^,'GRAMMA,I+ 

CALCULATE  ELEMENTS  OF  ROTATION  OPERATOR  ’  .  . 

.....  SLAT  =  SINF  (F.LAT/RAD) 

CLAT  =  COSF  (FLAT/RAD)  '  . . . -  -  -  - 

.  _  SLONG  =  S IN F  <  FLONG/R AD ) 

CLONG  ■■  COSF  (FLONG/RAO)  . * . 

_ TS.TART  «  TS  T ART+GBI  AS 

TSTOP  *  TSTOP+GBIAS  -  - - - - —  . . .  . 

~  READ  BINARY,  input  RECORDS 

READ  RUN  NUMBER  . .  . . -  .  . 

READ  TAPE  10,  IRUN 
IFINROCKT-IRUN  )45,  53,45 
45.  PRINT  50 

50  CALLA£)m/52H  WR0NG  INPUT  TAPE  WAS  M°UNTED  0N  B-5.  RUN  CANCELLED.) 

WRITE  RUN  NUMBER  . . . * . *  . 

WR  ITE  TAPE  9  » I  RUN 

READ  INPUT  OATA  . .  ‘  *  "  "  ’ 

READ-TAPE.  IQ, I  COUNT,  ( (OUT!  I.,  J),  1=1, 15),  J  =  KICDUNT) 

CHECK  STARTING  AND  STOPPING  TIMES  ~~  . . 

IE  (OUT.  1-U1J  -J5.TAR  T 1 62 , 65 , 50 

IF(OUT(l,l)-TSTOP)65,65,50OO  .  '  ~  ""  . 

IF  (OUT  ( 1 ,  ICO  UN  Tl-TSTAR.T)  55,65,65 

CONTINUE  . * 

DO  200  K=l, ICOUNT  . 

TL  =  OU T ( 1 , K )  -  GBIAS  . 

YM  =  QU  T ( 7 , K1 

XM  *  OUT (8  ,K )  . — . 
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62 

65 


OUT  ( 1 ,  K  )  =  OUT  (  l.KJ-BIAS 
CHECK  IF  TL  IS  IN  (TLA, TLB) 
IF(TL-TLA)72,71,70 

70  IF(TL-TLB)71,71,72 
TURN  SWITCH  2  ON 

71  ISW2  =  1 
T2  =  TL 
GO  TO  73 

c  .  TURN  SWITCH  2  OFF 

72  I SW2  =  2 
. T.2  =  TL 

CHECK  IF  PAYLOAD  ABOVE  LIMIT  OF  ATMOSPHERE 
.73  IF(TL-TA)74,77,77 
SWITCH  TWO 

74  GO  TO  (75, 76), I SW 2 

75  T 1  =  TLA 
GO  TO  BO 

76  T1  =  TL 
GO  TO  80 

77  T 1  =  TA 
COMPUTE  CA  VECTOR 
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C  DETERMINE  COORDINATES  OF  CONING  AXIS  VECTOR 
80  CALL  TRAJ(T1,XA,YA,ZA,XAD0T,YAD0T,ZAD0T,ERA,HA) 

C  A  < 1 )  =  XADOT 
CA ( 2  )  =  YADOT 

C  A  (  3  )  =  ZADOT . .  _ 

C  COMPUTE  VEHICLE  POSITION  ANO  VELOCITY  AT  GAMMA  SCANNER  TIME  T 
CALL  .  TRAJ!  T2,X,Y,Z  ,  XDOT ,  YDOT ,  Z  DOT ,  ER ,  H ) 

C  ROTATE  FROM  J.I.  TO  GEOCENTRIC  COORDINATES 

ZSUM  =  Z+RE 
ZP  =  ZSUM*CLAT-Y*SLAT 
XP  =  ZP*CLONG-X*SLONG 
YP  =  ZP*SLONG+X#CLONG 
ZP  =  Y*CLAT+ZSUM*SLAT 

C  COMPUTE  LATITUDE  AND  LONGITUDE  OF  PAYLOAD  AT  TIME,T 
ELAT  =  ATANF  (ZP/SQRTF  (  XP*X P+ YP* YP )  ) 

ELONG  =  ATAN1  (YP,XP) 

C  .  .  CALCULATE  ELEMENTS  OF  ROT  ATI  ON .  OPERATOR.  _ _  _ 

CONI  =  E  LON  G-F  LONG /RAO 
CON  2  =  ELAT-FLAT/RAD 
CON  3  =  SINF(ELAT) 

CON  4  =  COSFIELAT)  . . 

C  DETERMINE  EAST, NORTH  ANO  VERTICAL  MAGNETIC  FIELD  COMPONENTS 
CALL  MAGNEJIH,  ELAT,  ELONG,  .0M.,  BR ,  BTHET  A.fiPHI)  _ 

C  COMPUTE  FR  VECTOR 

C  ROTATE  _IN TO  J.I.  CUOROINATES  .FROM  GEOCENTRIC  COORDINATE  -SYSTEM  AND 
C  THEN  CALCULATE  THE  DIRECTION  COSINES  OF  THE  MAGNETIC  FIELD  VECTOR 

FR(1)  *  ltJPHI+CUNl*(C0N4*BR-C0N3*BTHETA)  l/BM  .. 

FR ( 2  )  =  (CON1*SLAT*BPH1+BTHETA+CON2*BR)/BM 
FR1.3I  =.  ( -CONlsCLATY  SPHI-C  0N2*E  THETA  +  BR)  /BM 
CALL  ATUDE 

C  DETERMINE  .DIRECTION  COSINES  OF  DETECTORS  IN  PAYLOAD  HIIH 
C  RESPECT  TU  EAST-NORTH-VERTICAL  SYSTEM 

DO  85  I  <=  1,3  ...  ....  . 

COSDXUI  =  A(  1,  1 )  *TRIG1  (  I  l+A(  2, 1  )*TRIG2(  I  >+A<  3,1  >*STHETA(  I  ) 
COSOY(I)  =  A( 1,2)*TRIG1( I  >+A( 2,2)*TRIG2I  I  >  +  A  <  3  ,2  )  *S  THE!  A  ( I ) 
COSDZ(I)  =  Al 1, 3)*TRIG1 (1)+A( 2,3)*TRIG2U )+A<  3,3 )#STHETAII) 

C  COMPUTE  DETECTOR  ATTITUUES 

DPHIIII  =  ATAN1  ( CO SOX ( I ) , COSDY ( I ) ) 

85  DTHE1  A  ( I  )  =  ATANF  (  COSDZ  ( I  ) /SQRTF  ( 1. -COSDZI  I  )  **2.) ) 

C  OUTPUT  ROUTINE 

C  GENERATE  THE  46-ITEM  OUTPUT  RECORDS 

C  VERTICAL  GAMMA  SCANNER ( 20  AND  90  DEGS)  IN  CUUNTS/10  MSECS. 

C  HORIZONTAL  GAMMA  SCANNER! 20  AND  90  DEGS) IN  COUNTS/10  MSECS 
DO  1086  1=2,5 

1086  OUT  (  I  ,  K I  =  0UT(I,K)*l.E-2 
C  POSITION  OF  PAYLOAD  IN  KILOMETERS 

OUT!  16, K)  =  X 
OUT (  17, K )  =  Y 
OUT!  18, K.)  =  Z 

C  ALTITUDE  OF  VEHICLE  IN  KILOMETERS 

OUT ( 19, K)  =  H 

C  THEORETICAL  FIELD  COMPONENTS  IN  J.I.  COORDINATES 

OUT ( 20 ,  K )  =  FR!  1)*bM 
OUT ( 2 1 ,  K  )  =  FR(2)*bM 
OUT ( 22, K)  =  FR(3)*bM 

C  DIRECTION  COSINES  OF  X,Y,Z  MAGNETOMETERS 

L  =  22 
00  86  1=1,3 
DO  86  J  =  l,  3 
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o  n  o  o 


L  =  L  +  l 

86  OUTfL.Kl  =  A ( I ,  J ) 

AZIMUTH  AND  ELEVATION 

VERTICAL  GAMMA  SCANNER  ( 20DEG)  IN  DEGS.,  OUT  I  32 ,K ) , OUT  I  35 , K ) 
HORIZONTAL  GAMMA  SCANNER  <  20DFG)  .  IN  DEGS  .OUT  ( 33 ,K  I , OUT  (  36  , K  > 
BETA  DETECTOR/PHOTOMETER  IN  DEGS.,  OUT  ( 34 ,  K) ,  OIJT  (  37 ,  K ) 

DO  87  1=1,3 
L  =  L+l 

OUT  (  L  ,  K  J  =  DP,-II<I  )  *R  AD 

87  OUT ( L+3 , K )  =  DTHETA ( M *RAD 

-  C  RANGE  (BURST  TO  PAYLOAO)  IN  METERS 

OUT ( 38 , K )  =  SORTF ( (X-XB )**2+( Y- YB>*#2+ ( Z-ZB 1**2 ) 

C  AZIMUTH  AND  ELEVATION  OF  PAYLOAD-BURST  POSITION  VECTOR  IN  DEGREES 
OUT  (  39  ,  K)  =  R  AD<‘AT  AN1  (X-XB,Y-YB) 

C  OUT ( 40 , K  )  =  RAD#ARS1NI (Z-ZB) /OUT(38,K) ) 

XX  =  Z-ZB 

OUT ( 40 , K  !  =  RAD«ATANF(XX/SQRTF(0UT(38,K)**2-XX**2) ) 

C  AZIMUTH  AND  ELEVATION  OF  ROCKET  AXIS  IN  DEGREES 
OUT (41 , K  )  =  RAD*ATAN1(A(3,1),A(3,2)) 

C  OUT  (  42 ,  K I  =  RAD#ARSIN(  A<  3,  3)  > 

OUT  (42 ,  K  I  =  RAD*ATANF(A(3,3)/SQRTF(1.-A(3,3)**2.)  ).. 

C  PHOTOMETER  FUNCTIONAL  VALUE  IN  WATTS/(CMS*»2)/STERAD1ANS 

OUT (43, K  )  =  0UT(12,K)*P0TMET  . . 

C  BETA  DETECTOR  FUNCTIONAL  VALUE  IN  MEV/CM#*2/SEC~ 

1  =  1 

NBVLTS  =  10.*0UT(10,K) 

IF  (NBVLTSm.94,90 
90  IF(NBVLTS-NBETA)96,92,92 
.92  I  =  NBETA  .. 

94  OUT ( 44 , K I  =  ENERGY ( I ) 

GO  TO  100 
96  B VOLT S  =  NBVLTS 

OUT  ( 44  ,  K  )  =  BNERGY  (  NBVLTS  1  +  (  BNERGY(  NBVLTS+  1 1-BNERGY  (NB-VLTS)  )  * 

1<  10.*0UT(  10,K)-BV0LTS) 

c  TOTAL  GAMMA  DETECTUR  FUNCTIONAL  VALUE  IN  MEV/  (  (CM*<‘2  ) *SEC  I ) 
1001=1 

IF (OUT <11,K) -VOLTS! 11)130,130,110 
110  DO  120  1=2  iNGAMMA 

IF (OUT ( 1 1, K) -VOLTS! I  1)140, 130, 120 

120  CONTINUE  . *  •  • 

I  =  NGAMMA 

130  OUT (45, K)  =  VOLTS! I ) 

GO  TO  200 

140  OUT ( 45 , K !  =  GRAMMA  1 1-1 ) +C OUT! 1 1 , K)-VOLTS< 1-1 ))7(  VOLTS  I i ) -VOLTS ( 1 
1-1 )  1  *  ( GRAMMA!  I )  -GRAMMAl  I  -1 )  ) 

C  ANGLE  BETWEEN  BETA  DETECTOR  AND  MEASURED  FIELD 
C  OUT  146  , K  )  =  RAD*ARCOSF< XM/( XM*#2+YM#*2+ZM4*2) )  . 

200  0UT(46,K)  =  RAD*ATAN1  (  SQR  TF  I  (  XM»*2+  YM4*2+  ZM*<>2  )  #*2  -XM**2 )  XM  ) 

C  END  OF  MAIN  LOOP 

C  WRITE  BINARY  OUTPUT  TAPE 

WRITE  TAPE  9 , 1  COUNT ,  ( ( OUT  (I,J),I  =  1,46),J  =  1,  IC  OUJVT  ) 

KKK  =  (  IC0UNT-IC0UNT/50*50  +  49)/50+IC0UNT/50 
00  240  K=  1,  KKK 
N  =  50*K 
M=N-49 

N  P  AGE  =  NPAGE  +  1 
I  =  NPAGE 

WRITE  OUTPUT  T APE  1 2 ,  2 10 , NR OCK T ,  I ,  (  OUT (  1 ,  J )  , OUT  ( 16  , J  ) , OUT  ( 17  ,  J )  ,0UT 
1 ( 18, J ) ,OUT! 19, J) , OUT ( 10, J) ,OUT ( 1 1 , J ) , OUT ( 12 , J 1 ,OUT (8 , J ) ,OUT( 7, J ) ,0 
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- ..  .  2UTJ  9  y  J  )_»0U  T  (  6 1  J  ) » OUT  (  14»  J ) r  OUT {  13. J),  OUT  ( 1 5  •  J  )  •  1«M  hi  \ 

210  i2TT  1  1H18X6HB00K  1 22X 53HOUf PUT  FROM  A T f  It UD EUFTETfi iTTSTTOKl — P Rfffil- — 

-  . J  AM  FOR  ROCKe  TI  .3 ,  18X4HP  AGE  I  5///2  5X19HPA  YLOAD  POS  IT  ION  IN16X  16HDETFC 

—  - rib ~r---!:r e  ta_tota L  _P HJTO-  6  X  8H  RE  A  D I N  GS  8  X3  H AGC  4  X 1 4Hn I F  F 

AGAIN  0/25X1 HX8X1HY8X1HZ9X1HH 11X11 HGAMMA  MET  ER4X2  HXM4X  "’HYM4XPW  7M - 

. — "  l  biTfiHiflsT  ‘  ^ — 

*ttm£q3o6£TA  0etect°K  ORIENTATI  0N7X21H20.  D£Gl*  _GAJ1MA  ^CAhlN£B /1 1X6HH 

_ .  1  "““rS*  ■ 

8J)/y F  15  S!FjriDS74r2«I?,?X^G,4X5H,0E®,8X5H,0^^HT“G'‘ 

rfl'Xrt  nnT^f’230’NR0CKT’  i»(out'1,j'>,qut(3,J^OUT(5^jV,OUT(2  ' 

~23QT  4^!?$HBQDK  322AS3HflUIPUT_EJMMl.ATTmlnF  i>PT»M,u,r,n., . 

5crmoR  RDEXETI3'  18X4HPAGEI5///24X21KGAMMA  SCANNER  OUT PUTS1 1X26 HDFT 
2ECT0R  FUNCTIONAL.  VALUES8X28HTHE0RET  ICAL  F LE.LD  XDMPONFNT  R/i  i  yZu^tI 
3ME6X9H90  DEGREE 5X  9H 20  0EGRE£7X4HBETA5X5HTSi^S^!f  CO^dIn 
- ^  S/22X5  HMQBJ  Z3X.4.HVE  R  T2X.5  HHORI Z  3  X4  B  VJE  RT14X5HGAMMA6X 

6 10^M5EC^9X15H^MEV/ CM4«  2/SFC^  PYi7w3HRR'3  '^^^^M^S^C^Toiri^HTcOUNTs"/  " 

240  CONTINUE _ 

C  READ  IN  NEXT  DATA  RECORD  '  - 

GO  TO  55.  _ 

C  WRITE  END-OF-FILE  ON  OUTPUT  TAPE  - 

5  000  END  -Fi LEG¬ 


END  FILE  12 


END  .FILE  .13 


CALL 

REWULD 

(9) 

CALL 

REWULD 

(  10) 

CALL 

REWULD 

1  12) 

CALL 

REWU-UD 

<  13  L 

CALL 

END 

EXIT 
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TABLE  3.1  INPUTS  FOR  TRAJ  SUBROUTINE 


Rocket 

a 

tQ  (sec) 
from 
launch 

zo 

(km) 

*0 

(km/sec) 

Po 

(km) 

i0 

(km/sec 

)  (km/sec2) 

8 

26°12 ' 

110 

137.740 

2.7770 

58.310 

1.6350 

5.46  x  10*4 

9 

23°30' 

110 

140.620 

2.9370 

40.710 

1.1686 

4.78  x  10*4 

15 

21° 

40 

30.520 

1.4190 

9.030 

0.4964 

0 

19 

135° 

30 

29  .830 

1.5790 

5.730 

0.2937 

0 

26 

113°30' 

35 

29.834 

1.5790 

5.729 

0.2937 

0 
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Z  (Vertical) 


X  (East) 


Figure  3.1  Definition  of  coordinate  system. 


¥ 


Figure  3.2  Angle  9  computation  (ATAN1 ). 


8 


Figure  3.3  Longitude  <p  computation  (ATAN1). 


z 


Figure  3.4  X-Y  computation  (ATAN1). 
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CHAPTER  4 
COMPUTER  RESULTS 

4.1  OUTPUT  FORMAT  AND  DEFINITIONS 

The  machine  computations  and  digitalized  telemetry  data,  as  output 
from  the  computer  program,  were  printed  from  tape  onto  three  copy  com¬ 
puter  printout  papers.  Three  data  books  for  each  rocket  were  required 
to  accommodate  all  of  the  data  items.  Samples  of  the  printout  format 
are  Included  in  Tables  4.1  through  4.3.  The  format  was  con¬ 

structed  with  title  headings  for  each  column  of  data  to  facilitate 
easy  data  analysis.  Time  in  H-time  is  printed  on  the  lead  column  in 
each  book  in  seconds,  to  the  nearest  tenth  of  a  millisecond.  The  data 
interval  is  approximately  11  milliseconds. 

The  telemetry  data  and  computed  data  from  Books  1  through  3  (Tables  4.1 
through  4.3),  reading  from  left  to  right,  are  interpreted  aB  follows: 

a.  Payload  position  in  Johnston  Island  coordinates:  X,  Y, 
and  Z  in  kilometers,  as  defined  in  the  earlier  discussions, 
are  the  east,  north,  and  vertical  displacements  in  a 
Cartesian  coordinate  system  based  on  the  geographical  loca¬ 
tion  of  Johnston  Island  (Table  4.1). 

b.  Altitude  in  kilometers  is  the  vertical  distance  to  the 
earth's  surface. 

c.  Detector  outputs:  beta,  total  gamma,  and  photometer  are 
the  unconverted  telemetry  data  froni  each  detector. 

d.  Magnetometer  readings  Xm,  Y^,  Z m  are  given  in  gauas  units 
as  taken  from  the  digitalized  telemetry  data. 

e.  AGC  is  the  automatic  gain  control  given  in  microvolts  for 
Rockets  8  and  9  and  decibels  below  a  milliwatt  for  all  other 
rockets.  Values  above  0.7  microvolt  or  110  dbm  indicate 
poor  data  quality. 
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f.  DIFF  is  the  difference  between  unity  and  the  sum  of  two 
squared  quantities.  Those  quantities  are  the  ratio  of  the 
smoothed,  gain,  and  bias  adjusted  magnetometer  value  to  the 
maximum  magnetometer  value,  for  x  and  y  magnetometers, 
respectively.  The  number  is  an  index  of  angular  error 
involved  in  the  attitude  angles.  As  error  increases,  the 
absolute  value  of  the  number  increases. 

g.  Gain  ratio  is  the  ratio  of  the  maximum  x-magnetometer  reading 
divided  by  the  maximum  y-magnetometer  reading  and  is  com¬ 
puted  over  a  100-data-point  group  of  digitalized  magnetometer 
values.  The  number  is  again  an  indication  of  angular  error 
as  it  differs  from  unity. 

h.  CODE  indicates  a  time  gap  in  the  digitalized  data  by 
bracketing  asterisks. 

i.  A,  B,  and  C  (Table  4.2)  are  classified  numbers  which  are  defined 
in  Reference  2. 

j.  Rocket  axis  orientation  azimuth  (degrees)  and  elevation 
(degrees)  are  defined  as  the  angle  taken  positively  from 
north  and  positively  from  horizontal,  respectively. 

k.  Beta  detector  azimuth  (degrees)  and  elevation  (degrees)  are 
the  pointing  direction  of  the  beta  detector,  with  azimuth 
taken  as  the  positive  angle  clockwise  from  north  to  180 
degrees  and  negatively  counterclockwise  to  -180  degrees. 

The  elevation  angle  is  positive  above  horizontal  and  negative 
below  horizontal. 

l.  Orientation  angle  to  the  earth's  field  is  given  in  degrees 
and  is  the  angle  between  the  beta  detector's  pointing 
direction  and  the  earth's  field  vector. 

m.  20-degree  gamma  scanner  orientation  gives  the  20-degree 

and  90-degree  horizontal  and  vertical  gamma  detector  pointing 
directions  in  azimuth  and  elevation  in  degrees.  The  title 
is  incorrect  and  should  properly  be  90-degree  gamma  scanner 
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orientation.  20-degree  gamma  scanner  orientation  is  ±180 
degrees  from  the  indicated  angles.  The  angle  convention  is 
the  same  as  that  of  the  beta  detector. 

n.  Gamma  scanner  outputs  are  the  telemetry  data  in  counts  per 
10  milliseconds  (Table  4.3). 

°*  Detector  functional  values  give  the  converted  (volt  to 
indicated  units)  value  for  the  beta  and  omnidirectional 
gamma  detector,  in  MeV  per  square  centimeter  per  second, 
and  the  photometer  outputs  in  watts  per  square  centimeter 
per  steradian. 

p.  Theoretical  field  components  in  Johnston  Island  coordinates 
are  the  outputs  from  the  MAGNET  subroutine  rotated  into 
Johnston  Island  coordinates.  FR(1)  is  the  eastern  component, 
FR(2)  is  the  northern  component,  and  FR(3)  is  the  vertical 
component . 

A. 2  INPUT  DATA  DESCRIPTION 

The  computer  printouts  for  each  rocket  include  an  Input  Data 
Sheet  which  includes  all  initial  conditions,  time  inputs,  start  and 
stop  times,  trajectory  information,  and  beta  detector  and  omnidirec¬ 
tional  gamma  detector  conversion  tables.  These  sheets  are  presented 
in  Tables  4.4  through  4.8,  as  taken  from  the  final  printouts.  Table  4.7  is  not 
final,  since  the  first  attempt  to  reduce  the  data  for  Rocket  19  was  unsuccess¬ 
ful;  the  table  represents  only  the  first  attempt. 
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TABLE  4.1  OUTPUT  FORMAT,  ROCKET  19,  BOOK  1 
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S  M« 


o  O  —  —  5  —  *  *  «  £  9  o  —  £  ^  *  *  *  *  2  »  o  —  X  *  *  *  «  a  *  o  —  *!  2  *  £  2  2  2  2  O  —  *  *  *  $  ®  «  ♦  ©  —  *!  S' 

a  3  k  9  ®*  V  <P  999  ©  O  O  O  O  O  O  O  U  O  —  —  —  —  —  —  —  —  —  *4  *4  *4  *.  *4  4*  4*  4*  4*  **  m  m  m  m  m  m  m  <*  —  *  *  ♦  ♦ 


8  ra 


r  n  k 

Ml 

5? 


—  4s..»-.#tf4,r««»»>c*?»co  —  ''.--♦♦♦'/>««4'*>k®*?oo-  —  -  *,  —  ♦  ♦./'♦««*>f>«D<e 
•  ••♦♦♦♦♦♦♦♦♦♦♦OOCOCOOOO-  —  —  —  —  —  —  —  —  —  AAAN  AN  NRN»  #  i*«l*  r  A  A 


*£’ 

O  S 

« 

o  • 


t  I  I  t  I  I  t  I  •  •  t  I  •  I  I  I  t  I  I  I  I  t  I  I  I  I  I  t  I  I  t  I  I  t  •  t  t  I  I  »  •  I  I  I  t  I  t  I  I  I 

»*.**•  *  —  -909—  *>*®o<*.*«o  —  *■•'©*.*>'■-**..#4'*  —  *  «  «  —  4*«»o*'*«o*‘***-o»'4.#*- 
o  —  —  —  —  *4*4*4*41*  —  »  —  .#***,r**,r,r,«*««<>*.4»*.*.««««-e*»*c>o©o— — —  —  *4*44*4* 


►  #*NO 
♦  —  *4  —  O 

1  9  1 

’  o  ’ 

-  ♦  r-  9  -  ' 

o  o  o  o  o  ■ 

—  c 
o  * 

Sg 

1  **«**«"•  ' 
'*■♦•#  4*  — 

* 

cr 

«■“  *  *  4*  —  <* 

*•  —  o  *  •  < 

ff  ff  F  f  *  * 

»*  «  *■  I 
®  «  «r 
I*  ♦  4*  ' 

«  «  •  ' 

<*  *.  o 
^  —  o 

«  ®  • 

9  i 

4T  1 

*■«#  —  —  o®« 

♦  —  —  o 
«*-«*> 
«  «  «  « 

«. 

«r 

« 

*» 

*4 

o 

*  *»  1 

<r  ■*  • 
—  O  ' 

♦  * 1 

a.  — 
*.  — 

*  • 
*  ♦ 

£ 

<*«*•» 

-  -  - 

4 

4*  4*  —  4* 

- 

— 

4*  4*  i 

4  4 

1*  4* 
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TABLE  4.3  OUTPUT  FORMAT,  ROCKET  19,  BOOK  3 


o  o  o  o  o  o  o 

A-UiAiA 

-Off 


0  0  0-30 
I  I  I  I  I 

4  C  N  4  O 

(>•  .O  >4  —  9 

rvi  —  o  9  r~- 

>4  4  4  W 


tfMAir  IT  4MT  IT  a 
9999999C 
NNNNNNNf 


>000000 

11(111 

)4N-ffOl" 

4  •>  4  N  -  C 

k-  r»  ^  k-  k>  k- 

«r  4>  tr>  ir 

9  9  9  9  9  9 

1NNNNNN 


o  o  o  o  o  o  c 

I  I  I  1  (  II 

—  —  rfN  O  <4  *  C 
OK  4  N»  o  a 
9  f*  *  *(  c'  k.  - 

MMNNNNf 

r-  r-  f-  c-  f 

a  9  9  9  9  9  e 
N  MM  INN  N  f 


—  ♦  «  o  «•>  o  o 
nOMTKff  K 

MN---00 

o  o  o  o  o  o  o 

9  9  9  9  9  9  9 


f*  K-  O  <4  K- 
>4  —  9  «  ** 
0  0  9-99 
0  0  9  9  9 


—  *  K-  O  >4 
-»lM«OK 
4C»®<DK 
9  9  9  9  9  9 


O  —  (/(  9  •*>  «f»  C 
»T  «N>  9  <C  >4  —  « 
K»  K-  <0  -O  -O  <0  w 
9  9  9  9  9  9  « 

«D  SC  SD  SC  «  ®  t 


oooooooooc 

AAA AAAAAi; 

ff 4^-sffnoaa 

SC  f-  4)IT«"N  — 


>irtrif>irirtr»u 
r  9  9  9  9  9  9  0 
aNNNNNNC 


)  m  o  C  NO  > 

'  tr  *  -4  ♦  >4  -sr  t> 

►  9  9  9  9  9  9  C 

;  SD  v  sr  sc  v  sc  a 


3  O  O  O  O  O  C 

LAA.UAJ 

94-9«4> 

t  r-  <  4  <"N  . 
►  9  9  9  9  9  fl 
D  <0  -C  O  O  -C  4 
IMIMA  ff>  If*  ff>  a 
h  9  9  9  9  9  0 
g  N  N  N  N  N  " 


►  Kg  O  9  *»>  ■C  < 
5  SO  iT>  fM  O  K-  a 
-  O  O  O  O  9  G 
C  9  9  9  9  *>  S 


99999999999999999 

OOOOOOOOOOOOOOOOO 


99999999999999999999999999*0 

OOOOOOOOO00000000000000000CC 


■  J  NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN 

jE  oooooooooooooooooooooooooooooooooooooooooooooooooo 

T  I  -  o  I  I  I  |  I  I  I  I  I  I  I  I  »  I  I  I  I  I  |  |  I  |  I  |  |  I  I  I  I  I  I  I  I  I  I  I  I  I  (  I  I  •  I  I  •  •  I  I  •  • 

u  ^  ■  .  |  ,,,  ,,,  ,,,  ,,,  ,,,  I,.  I,,  in  u.  ui  ui  u.  ui  ui  uj  u.  ui  iij  m  uj  uj  lij  iU  Ui  UJ  li*  UJ  UJ  Ui  U<  W  UJUiUlUUaUJUiUJttJUJ  (41  (4a  (L  (4J  (4J 

"5  Z  ooSoooooooooooooooooouooooooooooooooouoooooooooooo 

(  q  «  t/)  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOUOOOOOOOOOOOOOOO 

2^  U  IN  4'44444  4'  4444444444444444444444444444444444444444444 

(->  • 

£  &  >4  9  9  9  *  9  9  >4  ♦  9  *  ■»  •»  ♦  +  +  +  >4  <4  *  -4  •»  -4  -4  •»  -4  -4  >4  *4  -4  -4  -4  ♦  >4  >4  «4  >4  >4  4T  >4  >4  4T  «4  >4  >4  -4  >4  -4  >4  >4 

Oi  (J  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o 

°i  >  uiUjiiiUiiiJIUUJIuUJWU.U.iUaUJU'WUaUaWUdbUMUluU'UIUUU.UJU'UWWU.UaUWaKbttMUUMUB'UaUlUaiUUa 

Ul  (li  oooooooooooooooooooooooooooooooooooooooooooooooooo 

g  X  (^NNNNNNNNNNNNNNNNNNNNNNNNNNNNNIMNNNNNNNNNNNNNNNNNNN 

NCgfgNNNNNNNNNNNNNNNfgNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN 

►_  —  90<4*'SC94  —  IN*N9  —  9(N(N*ir(N®— (T9Ka»NSD>4^'0«"4iN©ir'>4«494'>4*N 

K  —  —  SD*^K*—*N«'a  —  O'O-O  —  —  ■O'ff—  -C-©  —  —  4'(ff'*'"-'©4'^CN  —  (/(►-■#< 

iA  UJ  >  oooooooooo-o  —  o-  ooooooooooooooo  —  ooo-ooooo  —  ooo  —  ooooo 


3  0  — 

0.  Ua  u 

(•on  ua  tr  «  —  tf>  — 

3  —  i/»  K(  ©  ©  4*  O 


°  £  8  *  OO'NOOOOOOOOINOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

«  z  o 

(11  — 

i  ^  ^  0*«««rN«**«K-K-**(K*'*  —  «  4  ■>  SD(r'ir®-#*«(r9»»(*K-«>*O-»9»-«^®*9S0»N*f2-» 

«  «  —  O— (NO  —  —  O  —  *•»  —  «  in  —  MTOiTO  —  —  ■C  -g>ra----in»i0Na>--NN0«0-»-®N-N4'C 

tri  ua  >  a  0000(NOOO*NOC>OOCOOO»NO(NOOOOOOOK<0000000000  —  0091000-»0000 

«  o 

«o  (-> 

S  Q  ig  *  N4("N04NN  WNNO«OK'^'lAON9-®(‘  N  9  ff1  N  4  O  4  K-09  N  N  >49  9  Ka 

«  h  N4CN44NIT  •NN®ff'®4'®«<®NO«94  N  —  4NO  ^9  ^  O  —  N  N  >4—  O  N 

^90  OOOOffwOOOOOOO  —  oc-oo  —  ooo  —  oooooo  —  —  —  OOOO  —  oooooooo  —  ooooo 
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»OTH€T  *  8.40E-08  W»TT/CB..Z/Srf K4DUN/V0H 


o 


3  UJU.UJlUUJUilLIUj 
-J 

tk  IPHIAnVi0KlA 


J  .  8  ....... 


3  UJliJUiWUUJUiUJUJ 
mi  (?NN-«0«JN0 

u.  o<o»N^#no 

iff  ooooooooo 

o  OHHNiNrt^f 


3  UiUjUJUJUMUUJIU 
U» 

•  •••««••• 


ooooooooo 

O<0N«^O«N 

•  ••••••«. 


g*  x  ODOOOPOO 

3  lu  lli  m  mui  m  »i|ll 

55  .t4  r  °°« 

tu  3  “*  **  ®  *0 

BU> 

£  ooooooooo 

H  nrviHNnfiinM 


H  ooooooooo 

3  UJUIlUUiUiUJliiUiUi 
J  f  ^«oarNi)D 
u. 

^  OOOOOOOOO 
H  N«<ro<ONSto 

§}  ou^NivnAtn 


3  W  Ui  Mi  IU  u  IU  |U  W  IU 
«  NNivai^awicf 

u.  op>p»n«ffl*«p>«» 

NNKMlf! 

*  WWUW wuugu 

5 
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TABLE  4.5  INPUT  DATA,  ROCKET  9 


UWUlilUJIUUJUl 

+  *~&'00'0-t0 


4-ir»«r»irin^-o^) 

oooooooo 


<  «M  f\  r»>  <  •#  IT> 


oooooooo 

UJUJIilULUUJUJUi 

U>«lft(rO<6<M<D 

OOOOOOOOO 

♦  o-c^w-ro-ora 

o  —  ■*  rg  inj  <«*■*■  tr* 


OOHftiNWim^  ia 


UJ  u.'  ID  UJ  IL'  U'  UJ  U  ’  UJ 

OBOOCMl 

cr,-*o«*rg^Mi^m 

OOOOOOOOO 

Mffi^-O'Orgc^'O 


UtfUJUJlUUJUJUJUJUJ 

^OnfflTM'IBO 


trs  ir  «  $  «  « 

K  O  ©  O  ©  O  © 

3  UJlUIUUliMU 
J  Nir»  004 

tb  tfVWN  ^  J 

n)hM^^no 


V)  rococo 

k  f*'  ir>  —  i»  g> 

^  oomnniao 


z  • 

O  X 

—  o 


<  K 

13 


UJ  Ui  UJ  UJ  UJ  UJ 
OIOMQ0K 

-<  40  W  «  _l  ^ 


OO-NNWO 


UJ  UJ  UJ  UJ  UJ  UJ 

BBlAHBCffi 

BOOinOBO' 


O  O  H  iJ  M  |*i  if 


UJ  UJ  UJ  UJ  UJ  UJ  UJ 
4B04#NB 
M  W  >©  •#  O*  <4> 

N  tf\  ff-  «-<  _«  in  4- 


O  O  ^  h  fvj  |*\  ^ 


UJ  UJ  UJ  UJ  UJ  UJ  UJ 
(M*-  O  •#  o»  rt  pi 
KBO'H'BIOO 


O  O  H  m  n  n  <r 


*AU<in4  B  4  J) 

o  o  o  o  o  o  o 


UJ  UJ  UJ  UJ  UJ  UJ  UJ 

®  o  t>  K  <?>  m  *4 

N»HNht# 


"ime 


»4NI«1 


2  ooooooo 

*“  «4f04(M(D«) 

^  OO-J^rgrgm 


j 
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CONTROL  PARENETERS 


D  U,  Ul  u UJ  Uj  IL'  UJ 
_i  <o  m  <nj  — 

u.  o  it  c  o  in  m  -t 


a  m  ^  , 


m  UJ  UJ  UJ  UJ  UJ  UJ  uj 


—  > 

D  O  — 
I  X  > 


-f  in  ui  ir\  ir  in  <c'«o 
OOOOOOOO 

H  Uj  UJ  U.  U.  Hi  Ul 1  UJ 
W-'f  0  0«\,'0 
MiMo-c-c-ri^r*- 


U  UJ  Ul  U  U  m  uj  1 1  il 
j-iiNiinoo^rvj, 


n  o  o  O 

<  uj  uj  Uj 

«  o  o  o 

I- 

O  IT  c- 
►-  O  eg  — 

uj  o  <r  n 

*  O  m  eg 

O  o  i*»  »r> 


3  UJ  UJ  UJ  Uj  Uj  UJ  UJ  UJ 

— J  NinNN  «o  «o 
u.  OH-roH»\r. 

je  -J  (\J  ■*  J3  CO  rJ  -J  o 


oo<-NNmn^m 


UJ  W  LU  UJ  Ul  U  W  UJ  i*  i 

r*ouoir»m«inoo 

«0>-»OBm  <■*  U3  rg 


oo-M\nm^in 


</i  oooooooc 

J  •«•••••• 

O  OO-HMOfifl 


O  IUUJUJUJUJUJUJUJ 

-J  + 

u-  g,o®o»n^«7'f*- 


UJUJUIUJUJUjUJUjUJ 

-n»NOrMO/0“ 

o«NCftr»ir>0*»-tf\o 


UJUJUJUJUJUJUJUJ 

00^«D«r^^n- 


8  oc:«HNrImin 
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TABLE  4.7  INPUT  DATA,  ROCKET  19 


*  to 

O  O  O  O 


UI  UJ  Ui  UJ 
—  N  I*  *\ 

•*  N  <*  O 


*a  eoooooo 

»-  ■ONI)'#  OlON 

J 

O  0-HNWW4 


UJ  Ui  Ul  UJ  Ui  U  UJ 

«  4  m  <#  o  o  o 

1*1  H  H  .4  N  N  W 

0N4OMNA 


^  0000000 

h  in  •*  r-  m  *  u\  -u 


o  O-J’JNNC'4 


S£ 


Ui  UJ 

5* 

o  z 
u  — 


on  OOOOOOO 
*-  ♦  O  *  N  «  *  O 


>r  *r»  m  it  o  o  co 
OOOOOOO 
K 

3  ui  ua  oj  uj  uj  uj  uj 

-J  BlflNtMnOO 

u.  n  >o  m  r- 


IflHinOrtK  H 


vn  oocoooo 

♦-  miMj'-Mnff' 

J  . ♦ 

a  ooHNNiflw 


*  mm  m  4  r>- 
OOOOOOO 
X 

3  UJ  Uj  uj  uj  uj  UJ  uj 

^  <tmOI*IHOO 

U.  S^OBNIfth 


m  Hin  U>  >^IA  h 


m  oouoooo 

H  N  B^U  lONOD 


o  oo-*nm*uo 


4AiniA«4f< 

OOOOOOO 

X 

3  UJ  ui  UJ  UJ  UJ  UJ  UJ 
-J  o  -<  m  oo  -o  m  an 
U-  NMNAOWN 

N  H  NlA  H  (*MA 


OOOOOOO 


O  OO-uftiw  ic\ 


K 

3 

aJ 

Ik 


mm.o*«r*K4-r»«j 

OOOOOOOOOO 


UJ  UJ  UJ  UJ  UI  UJ  U)  UJ  Uf  UJ 

K\>*«IJ(AO»«BN 


m  o>  «r  in  m  o  o  o  o  o  o 
k*  OA«»IA>*HA|hO 
J  #•••■••■••• 
O  OOOO-tMN'nmiPO 


OOOOOOOOOO 

UJUJUJUJUJUJUJUJUJUJ 

OUMn 

momwMnhoC'rj 

«J  o  o  o  o  oo  oo  O 
0n«»«04MBB 

oooo*<NNnn«o 


3  ui  ui  i^i  uj  in  in  in  i  in  in 

-J 

u-  -♦  o  iv»  _•  -  m  <o  r*  «m  — 

•M  4^N4BHN||UNO 

«/>  MflK'irOOOOOO 
I-  O  «N  m  »  <*l  CMCi  O 
"J  ••••••••••• 

O  OOOOH-N««4o 

> 


UJUJUJUJUJUJUJUJUJUJ 
off'OO'^o^uiMn 
O'  n^HnNrtftj#«o 


0000**Hi\imW40 


OOOOOOOOOOO 

X 

3  UJUJUJUJUJUJUJUJUJUJUJ 

-j  <n*Njr*O'0o*M«i-fuf'jr\» 

u.  «omn>*Bi*wn400 

0*MO'»-*^«0-*(Sj.r»-lrU 


ul  uMnuiirooooooo 

J  «•••••••••■ 

O  0  0  0'J"*rt(\|f<j(fl4ij\ 


•*  mmoo-or-r-r-r-co 
OOOOOOOOOOO 
K 

3  UJ  Ui  UJ  Ui  Ui  UJ  UJ  UJ  UJ  UJ  UJ 

_i  •rf-rnh'isima'-f.-amo 

U.  «Nm«kOHH»4ID 


m-uoo'uoj-o^rumr-  -o 


Ul  *  OOOOOOOOOO 
►“  ©^♦NC'flNB^OM 
«i  ••••••••••• 

O  O  OOOHHN  N«  4U» 
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UJtUUJUJliJlUUJUJ 

2SJSSSSS  _ 

•  ■"NI'IiOHhWO 


•#tririniri/N<04) 

^  oooooooo 

b  UlLUUJUJUtUJUJUl 

u. 

w)  ooooooooo 
O  o  «  h  n  rv  in  #  V  ^ 


3  UJUJkUUJLUUJUJUJ 
U.  ff'MIM-#  tr>  <r\ 

•  •»••••»  i 

^hNKI^Bi-No 

VI  ooooooooo 

*"  -rO«N*^o«N 
J  •♦•##•#•• 

O  1/1 


-n/NintTMAir-c* 

M  oooooooo 

3  UJUJUJUJiuuJUiUJ 
-j  ff-r  IMNNO 

u.  KhO(VO(M*Wvj 

Vi  ooooooooo 
^  ooJnlrlinwVin 


Stnini/>(/\ifi<0'Oo 

oooooooo 

Lii  U)  LU  U.  ui  L^i  g.i  ^ 
OOCO^NNOx 


Vi  ooooooooo 
►*  NBf  O-ONB^e 
!)  O  o  h  n  rj  m  m  ♦  i/? 


ooooooooo 

Ul  IJ  UJ  UMil  UJ  u  Ul  u 
*<OONOOHv 
m  ttr»  *•-««»««  * 


OOKXNHMH«« 


^  oooooooooo 

3  UJUJILIUJLiJUJUiUJUJUJ 
Uu 

o  —  N  i*.  <tc<Minco 


oo-hhn(v»i^^o 


4  H  •-  N  n!  ♦  f»'  H  r 


■i  *c  im  »  •*  a  o 

*  *4  +4  tsi  *\  +  4  O 


SsjisSisiSss 

*\  O'  HNi\rtf»n4lvii/io 


3  UJ  UJ  LDIUUJ.'UUJLDUJIU 

Vi  OOOOOOOOOO 
►“  u\  to  -*  ♦  n-  O  o  is;  ■>  * 

o  OO-HHfgrJinmVcj 


sss&sssssss 

3  U/UJUJtUOUIlUUittiUiUJ 

-i  vih>«Nhi 

«* 

v)  inmmmmnooooo 
§>  oo«-^NHNinfnV  V 


S<#irsiniAyMriiro«^ 

OOOOOOOOOO 

3  UiUJUJUIUJUJUJUJU^uJUJ 

■J 

IL  ■AMMN#l'ihiOOtOM 

£  90000000000 

8  o  o  **  »-  £  J  rv!  r>  w  V  ♦ 


3( 


